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ABSTRACT 


1 i i 


In February 1979 three Nike Tomahawk rockets were launched in Red Lake, 
Ontario, one previous to, and two during, the total eclipse of the sun, for 
the purpose of studying the atmosphere and it's interaction with solar 
radiation. On board were experiments from both the University of Illinois 
and the University of Bern, Switzerland. This paper describes the method and 
preliminary results of three of thes experiments that were used to measure 
solar radiation in the X-ray, Lyman-ot, and visible parts of the spectrum. 

The instrumentation designed for this investigation is described in detail, 
along with post-flight data processing techniques. The retrieved data have 
been processed to some extent to verify that a valid representation of the 
solar radiation has been obtained. 

The Lyman- a experiment yielded very good results, and preliminary data 
have been included in this report. The visible radiation experiment served 
as support for the other experiments and also functioned well. Due to a 
high level of energetic particles during the eclipse, the X-ray data were 
contaminated and are not presented. However, satellite observations 
indicate no appreciable lr/el of X-rays from the sun during the time of 
the eclipse. 
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1. INTRODUCTION 

This report describes experiments prepared for rockets launched in 
Canada, near Red Lake, during, and before, the total solar eclipse of 
26 February 1979. Of major interest is solar radiation in the UV and X-ray 
parts of the spectrum, both as they relate to the source (the sun) and to 
their interaction with the atmosphere. 

These experiments were part of an integrated payload also containing 
experiments for ion composition, electron density and energetic particles, 
with the broad objective of studying the chemical modelling of the lower 
ionosphere (which, here, because of the pre-determined rocket apogee, is 
defined to extend up to 135 km altitude). Thus the solar radiation 
experiments serve to give, to the modelling, the necessary input rates due 
to Lyman-a (121.6 nm) and hard X-rays fas represented by the wavelength 
range of 0.2 to 0.8 nm) . 

The importance of solar eclipses in the study of the ionosphere is 
further discussed in Chapter 2. 

The experiment to measure solar X-rays in the wavelength range 0.2 to 
0.8 nm is described in Chapter 3. Though this experiment performed well in 
the rocket flights, an unexpectedly high flux of energetic particles 
contaminated the data. Fortunately, other (satel ite) data show that the sun 
itself was not active and therefore that hard X-rays could not have been a 
significant ionization source. 

The Lyman-a experiment, described in Chapter 4, did perform well, and 
valuable data were obtained. In addition to the measurement of (energy) 
flux at this wavelength, used for calculation of ionization rates, these 
data are used in studying scattered Lyman-a, in observing the brightness of 
the solar chromosphere (as 3 function of radial distance), and in determining 
rocket aspect (when used in conjunction with the spin magnetometer). 

A sensor for visible radiation was also included in the payload, 
primarily as a check of the trajectory data and of the calculations of 
eclipse circumstances. This experiment, described in Chapter 5, also 
performed well. 

The integration of these experiments into the payloads for Nike Tomahawk 
rockets is described in Chapter 6. 
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The algorithms and progr,. .j for computer-processing the data from the 
Lyman-a experiments are giver, in Chapter 7. 

The eclipse circumstances of the rockets, calculated from the trajectory 
data, arc presented in Chapter 8, together with a preliminary evaluation of 
the performance of the experiments. 

Chapter 9 contains a summary of this report and some suggestions for 
additional data processing and analy. :s. 
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SOI AR ECLIPSES AND THE IONOSPHERE 

\ solar eclipse provides an opportunity to study the interaction of 
solar radiation with the ionosphere as well as the sun as a source of ionizing 
radiation. Since an eclipse provides a predictable obscuration of the solar 
disk, changes in the electrical composition of the ionosphere during an 
eclipse can be related to the ionizing radiation coming from the sun. 

Rocket investigations are suitable for this type of study because the solar 
radiation can be measured simultaneously with the electrical properties of 
interest at various levels of the atmosphere. This chapter will discuss 
work that has already been done in this area and results of previous rocket 
experiments . 

2.1 Electron Density 

During an eclipse, electron density measurements are made principally 
in the . • , and F regions of the ionosphere. In the F region (above about 
250 km) the eclipse effect is difficult to identify because of the generally 
less regular behavior and longer time constants associated with this region. 

The most regular of the ionospheric iayers, the .-layer (at about 
110 km), exhibits a clearly defined eclipse effect. The electron density 
is found to decrease gradually as the visible disk of the sun is obscured, 
but even in total eclipses of the maximum possible duration (approximately 
7 minutes) the electron density never decreases to less than one third of 
its normal value [L'mith at at., 1965]. 

Originally it was thought that ionizing radiation from the sun is 
entirely absent during totality and that the time constants of the processes 
by which the electrons are removed are such that the electron density does 
not fully decay to zero. . u has since been realized that the solar corona, 
which due to its size is not fully obscurrcd during a total eclipse, is a 
significant source of radiation so that even during totality there is a 
certain amount of ionizing radiation reaching the atmosphere causing 
reduced (but still present) electron concentrations. 

This theory lias been supported by previous eclipse rocket investiga- 
tions. Particularly during the eclipse of 7 March 1970 at Wallops Island, 
Virginia where Accardo et al. [1972] measured the residual flux of totality 
in the bands 0. 2-0.8, 0. 8-2.0, and 4. 4-6.0 nm to be 5, 7, and 16 percent of 
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the flux from the uneclipscd sun. Also Smith [1972] measured the Lyman-a 
radiation (1216 K) from the solar chromosphere and corona during the 
totality of an eclipse to be 0.15 percent of the flux from the unobscured sun. 

In the V region (nominally 60-90 km) electron densities have also been 
observed to decrease substantially below their normal value during an 
eclipse. [Smith et al., 1965; Mechtly et al., 1969; Mechtly et al., 1972]. 
This, as yet, has been not satisfactorily explained, partly due to the 
complex eher istry of the water-cluster ions in this region. This is the 
region of part'ular interest in the experiments for the eclipse of 
26 February 1979. 

2.2 Ionizing Rud\ition 

As the major source of ionization in the earth's atmosphere, the 
radiation generated by the sun covers a wide spectrum of energies including 
the more important X-rays and extreme ultra-violet radiations. The solar 
spectrum includes numerous emmission lines eminating from the chromospheric 
and coronal regions of the sun, as well as some continuum radiation. 

The present study is concerned with solar radiation which can penetrate 
to the lowest part of the ionosphere, the D region. The important ionizing 
radiation of the D region are the hard X-rays (0.25 to 0.8 nm) , the intense 
Lyman-a line emission (121.6 nm) , and the energetic charged particles that 
make up the galactic co ic rays [Whitten and Poppoff, 1971] . 

Measurements of solar radiation during an eclipse prove valuable in 
three ways. The first is that radiation measurements can provide information 
about the ’•adiation source; second, the concentration profiles of the 
constituents of the neutral atmosphere can be derived using absorption 
spectroscopy methods; and third, insight is gained into the mechanisms of 
radiation-ionospheric interactions by observing the variation in electron 
density caused by the obscuration of the solar radiation source. 

2.2.1 Solar radiation source. Previous eclipse flights have provided 
valuable information about the sun as a source if Lyman-a radiation. It has 
been confirmed that most of the solar Lyman-a radiation cminates from the 
chromosphere, which has been shown to extend approximately 4000 km above 
the surface of the photosphere (the visible sphere) [Weeks and Smith, 1971], 
The eclipse is thus of smaller magnitude for Lyman-a that it is in the 
visible spec. mi. This is illustrated in Figure 2.1 in which the crosses 
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give the Lyman- a flu.* as a percentage of the obscured sun as observed by 
Nike Apache 14.274 during the eclipse of 12 November 1966. The solid curves 
show the calculated percentage of the disk visible as derived from the 
trajectory of the payload for different disc sizes. The curve labeled a » 0 
refers to the size of k he solar disc as it appears in the visible spectrum, 
while the curves for a * 0.005, a = 0.01, and a = 0.02 refer to discs larger 
by 0.5, 1, and 2 percent, respectively. It is seen in Figure 2.1 that the 
chromosphere has a diameter approximately 0.5 percent larger than that of 
the visible disc. 

During the eclipse of 7 March 1970 further Lyman-a measurements were 
made with the aid of Nike Apache 14.437. It was found that the brightness 
of the chromospheric Lyman-a decreases exponentially over a range of 
0.005-0.03 solar radii from the photospheric limb, with a scale height of 
3835 ± 70 km [Smith, 1972]. Figure 2.2 shows this measured variation 
expressed relative to the brightness of the first 1 percent of the disc 
exposed, which was found to be approximately 14 percent greater than the 
average brightness of the entire Lyman-a disc indicating limb brightening 
of the chromosphere in Lyman-a. 

2.2.2 Absorption profiles. The measurement of the absorption profile 
of a certain wavelength of solar radiation is one method of determining 
concentration profiles of specific neutral constituents of the atmosphere. 
This technique requires a knowledge of the absorption cross section of the 
species in question over a frequency range inclusive of the pass-band of the 
radiation detector. An effective absorption cross section of the species 
can be derived by taking the weighted mean of the absorption cross section 
over the spectral width of the detector. When applying this technique one 
must be certain that the absorption of the species in question is indeed 
being measured and not some mixture of constituents with similar absorption 
properties. Measurements of Lyman-a ana 260 nm absorption profiles have 
provided data for molecular oxygen and ozone concentrations, respectively 
[Evans and Smith, 1975], Concentration profiles obtained during eclipse 
conditions provide information about the transport phenomena Oi the neutral 
species under investigation. 

In addition to deriving concentration profiles, absorption profiles 
are also useful in studying scattered Lyman-a radiation. 
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2.2.3 Ionization of the D-vegior.. As already has been noted the 
three main sources of ionisation of the 0- region (at about 80 km) are 
Lyman-a radiation, hard X-rays, and galactic cosmic rays. The very energetic 
galactic cosmic rays and solar X-rays ionize all atmospheric constituents, 
whereas the Lyman-a radiation can ionize only those constituents having 
ionization potentials lower than 10 eV. 

The relative importance of the three major sources of 0-region 
ionization was determined during the eclipse of 7 march 1970 by Aooardc et al . 
[1972]. It was found that at 80 km, during totality, the ionization of a 

trace constituent, nitric oxide, by Lyman-a produced 4 * 10 ion-pairs 

-3 -1 _"> 

cm s , the ionization by X-rays (0.2 to 0.8 run) produce 1.1 x io ion- 

-3 -1 

pairs cm s , and the ionization by galactic cosmic rays produced 
-2 . . -3 -1 

0.4 x io ion-pairs cm s , thus giving a total, for these three sources 

-2 -3 -1 

5.5 x io ion-pairs cm s . These data indicate that ionization of NO 

by Lyman-a radiation is the major source of ionization at the altitudes of 
the 0-region. This is because the intense Lyman-a line of the solar spectrum 
aligns with a transmission window in the absorption cross section of 
molecular oxygen (a major constituent of the atmosphere) allowing penetration 
of Lyman-a into the O-region, whereas the other ionizing radiations are 
absorbed at higher altitudes. Therefore, Lyman-a radiation is of major 
interest for interpretation of the changes in electrical properties of the 
O-region that occur during an eclipse. 
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3. SOLAR X-RAYS (0.2 to 0.8 nm) 

Each of the three payloads carried a Geiger counter sensitive to 
radiation in the range of 0.2 to 0.8 ran (1.5 to 6.2 keV) . The operation 
of the Geiger counter and associated counter circuit will be discussed, 
followed by the calibration procedures and calibration records for the 
X-ray experiment. 

Unit optical depths are 72 and 95 km for the wavelengths 0.2 and 0.8 nm, 
respectively. The band of X-rays selected for investigation is represent- 
ative of these unit optical depths and is a source of ionization in the 
D-region. The band of X-rays is convenient for measurement by a Geiger 
counter with a thin window of beryllium. 

Although this band of X-rays is relatively broad, the actual spectral 
energy distribution can be recovered by assuming a black body temperature 
for the source. This can be obtained from the absorption profile of the 
radiation, although not during the eclipse itself. The method is discussed 
in Aaaardo et at. [1972]. 

3.1 The Geiger Counter and Bias Cirau.it 

A geiger counter is an extremely sensitive device used for the detection 
of ionizing radiation, either corpuscular or electromagnetic. The high 
sensitivity of the Geiger counter is due to the fact that an individual 
photon or energetic particle can trigger a discharge from the Geiger 
counter tube. Thus any types of particles which can release charge within 
a Geiger counter can be detected, including X-rays which produce ionization 
by secondary processes [Price, 1964], 

Geiger counters generally consist of a chamber containing a conducting 
cylindrical surface which surrounds a wire along its axis of symmetry. 

A cylindrical geometry is commonly used because it facilitates parameters 
that are easily controlled to meet design requirements. The chamber is 
filled with an inert gas, usually helium, argon, or neon. A :mall 
percentage of a polyatomic gas is also added for quenching purposes. 

Quenching is the process of terminating an electrical discharge within the 
tube, which will be explained later in this section. Any polyatomic gas 
molecule will produce the quenching action although each polyatomic gas 
has certain advantages and disadvantages which must be considered when 
choosing a gas fill for a Geiger counter [Korff, 1947]. 


A*- 
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The Geiger counters used in the present experiment arc contained in a 
stainless steel cylindrical chamber which also serves as the outer conductor. 
The chamber has a wall thickness of 0.5 mm, radius of 12.5 mm, length of 
76 mm, a central wire of radius 0.005 mm, and is filled with neon at 0.5 Atm 
with 1% isobutane added as the quenching agent. 

For operation of the Geiger counter a large potential is applied between 
the outer conductor which is referred to as the cathode and the central wire 
which is called the anode or, sometimes, the collector electrode. This 
potential sets up an electric field defined by the geometry of the electrodes, 
with magnitude given by (3.1), 

E £ 7 T-^r (3.1) 

r l°g l0 (V r 2 ) 

where E is the electric field strength, V is the applied voltage between the 
cathode and the anode, r is the point of observation, and r and r* are the 
radii of the wire and the cylinder respectively. Using equation (3.1), the 
aforementioned dimensions of the Geiger tube, and an applied voltage of 
950 V (the voltage used in the experiment), the electric field intensity is 
found to vary monotonical ly from 79.2 kV cm at the wire surface to 317 kVcnf* 
at the cylinder wall. At only 1 mm away from the wire surface the field 
strength drops to 3.77 kV cm * indicating a rapid increase in the electric 
field when approaching the central wire. Keeping this in mind, an explanation 
of the discharge operation of the tube will follow. 

Assume that an ionizing event has taken place within the volume of the 
Geiger tube, i.e., one or more charged particles, or photons, or neutral 
particles have by some process produced ionization within the electric field 
inside the Geiger tube. Once this ionizing event has taken place the 
electric field will cause the positive ions to move toward the cylinder 
wall and the negative ions and the electrons to move toward the central 
wire. As the electrons approach the wire, they enter the region of rapidly 
increasing electric field and are therefore accelerated to high velocities. 
These electrons will cause the release of more electrons upon collision 
with any atoms in their path. The additional electrons released along with 
the original electrons will continue moving toward the wire causing even 
more electrons to be released along the way giving an avalanche effect. 
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Since the charge collected on the wire depends on the final number of 
electrons arriving there, a response due to one single electron is greatly 
amplified in this manner. 

While this process is taking place the positive ions, moving much slower 
due to their greater mass, drift toward the cylinder wall. Ideally, when the 
positive ions reach the cylinder wall they are neutralized terminating the 
discharge process and the counter tube is ready to detect another ionization 
event. In reality the process of neutralizing the positive ions may result 
in the release of electrons which will move inward and initiate new 
avalanches and cause the discharge process to repeat itself, resulting in 
a multiple discharge or a continuous, self-sustaining discharge. To prevent 
the undesirable release of electrons upon neutralization of the positive ions 
is the function of the quenching gas. The quenching gas will inhibit the 
release of electrons thus removing the conditions for reinitiation of the 
discharge process. 

The Geiger tube does not immediately recover from the discharge, in that 
the positive ions form a nearly cylindrical sheath around the central wire 
which disrupts the electric field and prevents the initiation of new 
avalanches near the anode. This situation continues until the ion sheath 
has migrated out of the high electric field region near the wire. The time 
interval during which the tube is completely insensitive to incoming radia- 
tion is called the "dead time", followed by the "recovery time" which is the 
time required for the pulse size to regain its original amplitude. Another 
time interval that is commonly defined is the "resolving time" which is the 
time after a count for which the entire measuring system is insensitive to 
incoming pulses. The resolving time is always greater or equal to the dead 
time. 

Figure 3.1 illustrates these time intervals for a Geiger counter 
operating at a high counting rate.* The resolving time of a Geiger counter 
presents an upper bound on the counting rate, but if the resolving time, t, 
is known, a corrected counting rate N can be approximated for any measured 
rate N by using (3.2) [O'Kelley, 1962]. 

*To actually "see" an output resembling Figure 3.1 would require a 
storage oscilliscope because each pulse is triggered at a different time. 
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Figure 3.1 Illustration of pulse shapes in a typical Geiger 
tube operating at a high counting rate. The dead 
time and recovery times are determined by the 
Geiger tube characteristics, but the resolving 
time depends on the electronic recording system 
[O’Kelley, 1962]. 
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Since the Geiger counters for this experiment are used to detect X-rays 
which have a relatively high photoelectric absorption cross section, a thin 
beryllium window is mounted in the side of the cylinder to admit the X-rays 
into the Geiger tube. Most of X-rays incident upon the window will enter 
the tube because of the relatively low photoelectric absorption cross section 
of beryllium and the thinness of the window. The window and gas fill 
combination determines the band of wavelengths for which the counter is 
sensitive. The efficiency of the tube in the low energy ~egion decreases 
because of absorption in the window and in the high energy region because 
of the decrease in the photoelectric absorption cross section of the gas 
fill. With the present beryllium window and neon gas fill, the counter is 
theoretically sensitive to X-rays in the range of 1.5 to 6.2 keV. However, 
during the eclipse it became evident the the counter is also sensitive to 
electrons with energies greater than 40 keV. 

3.2 Electronic Circuits 

The block diagram in Figure 3.2 shows the structure of the X-ray 
electronics. The resolving time, as described in the previous section, is 
determined by the duration of a positive going pulse from the output of the 
monostable multivibrator during which the circuit will not respond to any 
signal coming from the X-ray tube. The dead time of the detector is 
approximately 130 microseconds while the resolving time of the circuit is 
set to approximately 190 microseconds giving an extra margin to allow for 
any changes in the performance of the counter during the flight as well as 
variations in operating characteristics from one counter to another. 

The 950 V bias necessary for the operation of the counter in the Geiger 

region is derived from the +30 V power bus using a 3 watt regulated DC-to-DC 

converter module manufactured by Technetic, Inc. Figure 3.3 gives the 

details of the biasing circuit. The series combination of R. and C. serves 

4 4 

as a smoothing filter for the DC-to-DC converter output. The signal from 
the detector tube will have a 950 V DC component which is removed from the 
active differentiator input by the capacitor C^. The resistance isolates 
the power supply from the central wire, thus allowing the anode voltage to 
drop upon the collection of electrons following a discharge in the Geiger 
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counter. R-, also serves as a protective device in that the 1 Mfi resistor 
will limit the current to 1 mA if the counter should go into continuous 
discharge or if the high voltage portions of the circuit should be 
accidential ly short circuited during testing. Similarly Rj will provide 
current limiting if the capacitor C. fails in such a way that it acts as a 
short circuit. Neglecting because it is relatively large compared with 
R, and R^ and neglecting because it is relatively large compared with C, 
and the input capacitance C,, the time constant of the bias circuit is 
given by, 


T fl * RC - (R 2 ♦ R 4 )(R 2 ♦ R,) 


where the values for R 2 , R 4> C,, and C 3 are taken from Figure 3.3. Therefore 
the response time of the bias circuit is less than the 190 ms resolving time 
of the counter circuit and will not slow down the count rate. 


The sensitivity of the active differentiator is such that pulses of 
amplitude greater than the threshold of the pulse height discriminator will 
be produced for two consecutive particles coming into the detector not 
closer than one system dead-time apart. Thus a photon or an energetic 
particle entering a not fully recovered Geiger counter will produce a smaller 
than normal output pulse, as shown in Figure 3.1, but will be detected by 
the edge sensing characteristics of the differentiation circuit and will 


also trigger the monostable multivibrator via the pulse height discriminator 
given that the particle enters the tube at least 190 ms Jater than the 
previous photon or energetic particle. 

The pulse height discriminator performs a twofold function in that it 
provides the proper logic levels for the multivibrator input and also 
provides some degree of noise immunity by only reacting to inputs greater 
or equal to a predetermined threshold voltage. Figure 3.4 gives the circuit 
details for the pulse-height discriminator. The threshold voltage V T is 
selected by choosing the proper value for R 3 in the following equation. 


V 


T 


- 22200 



0.2 
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To transmit a series of pulses such as those coming from the multi- 
vibrator output would require a very large dynamic range. In order to 
present this information on the FM/FM telemetry system of the payload the 
following scheme is used, which decreases the dynamic range requirements 
of the telemetry c' ^nnel while retaining the information density of the 
detector output. 

As mentioned previously the monostable multivibrator will produce an 
output pulse of approximately 190 ws when a photon or energetic particle 
enters the counter. These pulses are fed into a 4 bit binary counter with 
output to a 4 bit digital- to-analog converter which will produce a lower 
bandwidth sig.ial consisting of a series of 16 step staircases as the counter 
cycles through its 16 different outputs. The signal is then amplified by 
the buffer amplifier to translate the staircase to a 0 to 5 V level 
compatiable with the telemetry channel. Therefore, to determine an average 
count rate over a certain time interval, one must determine the number of 
"steps" in that interval. 

3.3 Testing and Calibration 

Table 3.1 of this section gives a summary of the testing and calibration 
procedures for the X-ray electronics and are followed by Figures 3.5 and 3.6 
which give the circuit details and the component layout, respectively. These 
figures should be referenced while reviewing the testing and calibration 
procedure. Table 3.2 gives the calibration record for the Geiger counters 
flown in Nike Tomahawks 18.1020, 18.1021, and 18.1022. 

Figure 3.7 is a photograph of a fully assembled X-ray deck. The Geiger 
Tube can be seen in the lower left-hand corner of the picture. 
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Table 3.1 Procedure for the testing and calibration of the X-ray electronics. 

1. Adjust the DC to DC converter to give an output of 950 V. 

2 . Place a radioactive source near the window of the counter tube and 

observe the tube output at test point Tj to verify proper operation of 

the tube. The pulses should be like the ones shown in Figure 3.1 with 
a maximum magnitude of approximately 1 V. 

3. While the counter is counting select the capacitance C 9 that will just 
stop the damped oscillations at test point T, . Measure and record the 
magnitude of the pulses. The magnitude should be approximately 7 V and 
to achieve r his it is acceptable to allow a few cycles of damped 
oscillations if necessary as long as they damp out in less than one 
resolving time. 

4. Isolate the input to the pulse height discriminator and apply a sine 

wave to determine the threshold voltage. While observing the output of 

the monostable multivibrator slowly increase the sine wave amplitude 
until the monostable is triggered. Measure the record and the thres- 
hold voltage which would be the amplitude that just triggers the 
monostable multivibrator. 

5. Using the sine wave input to the pulse height discriminator as in step 4 

adjust the amplitude to a value greater than the threshold voltage and 

monitor the output of the multivibrator. Adjust R, to give a pulse 

o 

duration of approximately 190 ps. 

6. With input to pulse height disciminator as in step 5 adjust R g to produce 
a 0 to 5 V staircase at test point T g . 

7. Reconnect differentiator output to the pulse height discriminator, place 
test source near the Geiger tube window, and observe output at test 
point T g to insure overall proper operation of the Geiger counter. 




C? SELECT 
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1'igure 3.6 Continued. 







A photograph of the X-ray deck for the payload of Nike 
Tomahawk 18.1021. The larger aluminum can houses the X 
ray electronics and the smaller one contains the high- 
voltage bias circuit. The Geiger counter and aperture- 
defining plate are at the lower left. This deck also 
carries the solar sensor, shown at the upper right. 
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4. LYMAN-a RADIATION 

F.ach payload carries two photometers* sensitive to a narrow spectral 
band of solar radiation. Mien flux of the proper wavelength (Lyman-a) is 
incident upon one of these devices, a small current proportional to the 
intensity of the radiation is produced, which is fed into an electrometer 
(current-to-voltage converter). The following sections describe the operation 
of the photoionization chamber, the associated electronics, and the use of 
the two photometers in determining rocket aspect and Lyman-a flux. 

4 . 1 Fhotoionization Chambers 

A photoionization chamber is a device used to measure flux intensity 
within a specific spectral band, consisting of a chamber of gas, a window to 
admit flux, and an applied voltage across two electrodes in contact with the 
gas. A cross section of the photometers used is shown in Figure 4.1, where 
a voltage would be applied between the pin electrode and the outer shell. 

A photon of sufficient energy will pass through the window into the 
chamber, collide with a gas molecule, causing the release of an electron, 
thus ionizing the molecule. The ionized gas molecule is then attracted to 
the negative electrode, which is the center pin, and the free electron drifts 
to the positive electrode, which is the outer shell, producing a current flow 
in the photometer. The resulting current is proportional to the number of 
photoionized molecules, which is proportional to the number of incoming 
photons, which gives a linear relationship between current and flux intensity. 
Therefore, when calibrated, a current measurement gives the intensity of 
solar radiation entering the chamber. 

The upper and lower limits of the spectral band of a detector are 
determined by the ionization potential of the gas and the filter character- 
istics of the window, respectively. Nitric oxide, the gas fill used in the 
photometers, has an ionization potential of 9.25 eV, which implies that an 
incoming photon must have at least 9.25 eV of energy to photoionize a gas 
molecule. This gives an upper wavelength limit of 135 nm. The window 
material, magnesium flouride, will pass wavelengths of 112 nm or greater 
giving a bandpass of 112 to 135 nm for the N0-MgF o combination. 




*Thc names photometer, photoionization chamber, and photon detector will he 
used interchangeably. 


Figure 4.1 Cross section of ultraviolet radiation ion chamber f Evans and Smith, 1975]. 



2 "/ 


The NO-MgF., detector was chosen because of its narrow bandwidth around 
the Lyman-a line (121.6 nm) , giving discrimination against unwanted wave- 
lengths and for the desirable non-hygroscopic qualities of the MgF 0 window. 

4.2 Electronic Circuits 

The instrumentation that will be referred to as the Lyman-a electronics 
consists of one photo-ionization chamber, two stages of amplification, and 
an in-flight calibration circuit. A block diagram of the Lyman-a electronics 
is given in Figure 4.2. As indicated in the figure, the photometers have a 
nominal quantum efficiency of 60 percent at Lyrrian-a. This quantum efficiency 
indicates that 60 percent of the incoming Lyman-a photons will ionize gas 
molecules within the photo-ionization chamber. Since each ionization of a 
gas molecule results in a released electron, the quantum efficiency can be 
used to derive the sensitivity of the photometer. 

In response to an incoming flux of 1 photon s' 1 , a photometer with a 
quantum efficiency of 60 percent would produce the current given by, 


i _19 

I a = (1 photon s )(0.6 electron/photon) (1 . 60 x 10 C/electron) 

,--19 _ -1 ,.-10 . 

=10 Cs =10 nA 


Therefore, this photometer would have a sensitivity of Sp = 10 ^ nA/ (photon *s) 
The quantum efficiencies (as measured by the manufacturer! and corresponding 
photometer sensitivities for the six photo-ionization chambers of the present 
experiment, are given in Table 4.1. 

The current from the photometer is fed into the first stage of 
amplification which is provided by the electrometer. The electrometer is a 
specialized operational amplifier designed for the measurement of extremely 
small currents. It is configured as a linear currcnt-to-voltage converter 
with nominal sensitivities of 0.25 V/nA and 2.5 V/nA for the pre-eclipse and 
eclipse flights, respectively. 

The second stage of amplification, the buffer amplifier, provides 
isolation of the electrometer and also supplies current drive to the telemetry 
channel. The buffer amplifier is configured for unity gain for the pre- 
cclipsc flight and a gain of 10 for the flights during the eclipse. Referring 
to Figure 4.2, the nominal overall sensitivity for the pre-eclipse electronics 
is. 


PRE-ECLIPSE : 


G, =0.25 V/nA 
g 2 = 1.0 V/V 
ECLIPSE: G, =2.5 V/nA 

g 2 = 10 V/V 


+ 30 V 


LYMAN - a 
RADIATION 



TO TELEMETRY 
CHANNEL 


PHOTOMETER 


ELECTROMETER BUFFER 
AMPLIFIER AMPLIFIER 


Figure 4.2 Block diagram of Lvman-a electronics. 






Table < 

4.1 Sensitivities of Lyman 
(U = uncollimated and i 

-a instrumentation 
C * collimated). 


Q.E./100 

S x io 21 
P 

*FS * ‘O 3 

S x 10 
u 

18.1020 

U 

0.591 

0.946 

19.4 

0.244 

18.1020 

C 

0.621 

0.994 

20.2 

0.246 

18.1021 

U 

0.536 

0.858 

0.208 

20.6 

18. 1021 

C 

0.678 

1.08 

0.208 

26.0 

18.1022 

U 

0.610 

0.976 

0.205 

23.8 

18. 1022 

C 

0.640 

1.02 

0.188 

27.1 


quantum efficiency of photometer 
as specified by the manufacturer 

sensitivity of photometer 

(Q.E./100 e/photon) (1 . 60 x 10 coul/e) f 1 photon s' 

(1 photon s" 1 ) 

(Q.F../100) (1.60 x 10" 21 ) A/photon • s -1 


overall sensitivity of Lyman-a electronics (including detector) 


s A | f s Volts ) 

> P nhoton • / I ^ ) 


1 photon • s~* 1 
5 S 

p Volts 

*FS photon • s~ ^ 


full scale input current 




30 


* [10 10 nA/ (photon*s)] (0.25 V/nA)(l V/V) 

® 2.5 x 10 ** V/(photon*s) 

Also, derived in the same manner, the nominal overall sensitivity for the 

-9 

eclipse flight is 2.5 x 10 ' V/ (photon *s ) , which is 100 times greater to 
detect the reduced Lyman-a intensity of ccliprr conditions. Since the 
circuitry is designed to provide a 5 V maximum signal at the output of 
the buffer amplifier, the instrumentation has the capability to measure 
Lyman-a intensity up to 2 x 10^ photons s”* for the pre-eclipse payload 
and 2 x 10 photons s for the eclipse payload. Actual overall sensitivities 
and full scale currents for the six Lyman-a detectors can be found in Table 
4.1. 

An in-flight calibration circuit provides a reference current to the 

electrometer producing a 5 V buffer amplifier output from 40 km to 70 km 

using altitude-sensitive switches (baroswitches) . 

A detailed circuit diagram of the Lyman-a electronic^ and associated 

circuitry is shown in Figure 4.3. Determination of the component values 

marked "select" is discussed in Section 4.4. The 89.6 V bias voltage for 

the photometer is provided by four 22.4 V batteries (Mallory RM 412) 

connected in series. The Keithlcy 302 electrometer is chosen for its high 

12 

input resistance of 10 ohms and desirable stability characteristics. 
Manufacturer's specifications for the Keithley 302 are given in Table 4.2. 
Diode protection is used on the input of the 741 operational amplifier as 
a safeguard against any large transients which may occur during the flight, 
and at the output to prevent a negative signal into the telemetry channel . 

The MC7814/MC791 5 voltage regulators provide ±15 V using the ±30 V busses 
as power sources. 

4 . 3 Frequency Response 

The total Lyman-a waveform delivered to the telemetry channel is 
dependent upon the angular sensitivity of the detector, the spin rate of the 
payload, as well as the frequency response of the circuits. As already 
mentioned in the introduction, each payload carries two photometers used 
together for determining rocket aspect and for measuring Lyman-a flux. The 
details of the arrangement will be given later in this chapter but, for now, 


MC78I5 
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Table 4.2 Specifications for the Keith ley Model 302, 


Specifications measured at 25*C 

DC VOLTAGE GAIN. OPEN LOOP: 
Unloaded (Minimum) 

Full Load (Minimum) 

INPUT CHARACTERISTICS: 

Common Mode 
Resistance (Minimum) 

Shunt Capacitance (Maximum) 
Rejection (Minimum) 

Voltage Limit 
Between Inputs 
Resistance (Minimum) 

Shunt Capacitance (Maximum) 
Overload Limit 

Current Stability 

Offset 

vs. Time (worst case) 

vs. Temperature (worst case) 

vs. Supply (worst case) 

Voltage Stability 
Offset 

vs. Time (worst case) 

vs. Temperature (worst case) 
vs Supply (worst case) 

Current Noise 

0.1-10 cps (Maximum peak-to-peak 
Voltage Noise 
0.110 cps (Maximum rms) 

10 cps— 500 kc (Maximum rms) 

FREQUENCY: 

Gain Bandwidth Product (Minimum) 
Slewing Ra*e (Minimum) 

Rolloff (Nominal) 

OUTPUT: 

Amplifier 
Reference Voltage 
Regulation 

OPERATING TEMPERATURE: 
CONNECTORS: 

Input 
All Other 

POWER REQUIREMENTS: 

Voltage (positive and negative; 


MODEL 302 

12.000 

10 000 


10 ,? ohms 

5 picofarads 

±100 volts continuous 
±400 volts momentary 

10 14 ampere 
10' 1S ampere/24 hours 
Doubles every 5*C 
5xl0 15 ampere /% 

Adjustable to zero 
2 millivolts/week after 
1-hour warm-up 

150 microvolts/'C 
1 millivolt/% 

5xl0‘ 15 ampere 

10 microvolts 
100 microvolts 

150 kc 

0 1 volt/microsecond 

6 db/octave 

±10 volts (gi 5 milliamperes 


-0*C to +50*C 

Teflon-insulated feed through 
#20 pins 0.2* long, 0.2* goo 

15 volts regulated ±0.1% 


+Current (Plus output current) 
—Current (Plus output current) 
DIMENSIONS: WEIGHT: 


5 milliamperes 
5 milliamperes 
1* high x 1* wide x 
’/i6* deep; V5 ounce 
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we need only note that one detector is collimated and the other is not. 

The angular variation of the uncollimatcd detector is expected to approach 
a (cosine) " function while the collimated detector is expected to have a 
nearly triangular variation with angle. Figure 4.4 shows idealized waveforms 
from a collimated and an uncollimated detector, respectively. Hue to the 
spinning motion of the payload, these waveforms are periodic with a period 
7=1/ (spin rate) . 

Since an accurate measurement of the peak amplitude of the signals is 

necessary in the final analysis, the frequency response of the Lyman-a 

electronics must be considered to determine if any reduction of the signal 

amplitude will be introduced by the Lyman-a circuit. Referring to Figure 4.3, 

consider the Lyman-a electronics with the frequency compensation capacitor, 

C removed. The circuit is essentially two cascaded operational 
/ 4 1 

amplifiers. A first-order approximation to the transfer function for an 
op-amp with feedback will be derived. Consider the op-amp configuration 
in Figure 4.5 where and are resistors and is the effective 

capacitance in the feedback path introduced by the physical layout and the 
op-amp. Summing current at node a gives. 


v .-v 
a 


V -V 
o a 


fi l H C cff 

sltagc at 
following equation is also true, 


where v is the voltage at node a. 

a 


(4.1) 


Using an ideal op-amp approximation, the 



(4.2) 


Combining equations (4.1) and (4.2) gives the following transfer function, 


v 

0 

V . 

1 


*1 II ; off 


*1 ( rr ^V'off) 


(4.3) 


which has a single-pole roll-off at = 1 Since the Keithley 302 

inherently has a much lower gain-bandwidth-product than the 741C and 
considering their relative gains, it will be assumed that the frequency 
response of the entire circuit is that of the Keithley 302 (multiplied by a 
factor of 10 for the eclipse electronics) and, furthermore, it will be 
assumed that the response of the Keithley 302 is characterized by the 


A 
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transfer function in (4.3) with R. replaced by R D [ see Figure 4.3). 

1 t 

It will now be shown that by applying a step input to the l.yman-a 
electronics 3nd measuring the rise time of the resulting output, a value for 
C^. of the electrometer can be obtained thus characterising the frequency 
response of the circuit. The step response of the circuit is, 

V q = i’ ; . [1 - exp(-t/fl f C eff )], t > 0 (4.4) 

where v_. is the amplitude of the step input and ^ is the resulting output. 

Since rise time is defined as tic time it takes for v to reach 0.9 v., the 

o V 

following equation must hold, 


0.9 = [(1 - exp(-t /? //? f C eff )] 


(4.5) 


where t 


R 


for C 


off’ 


denotes rise time. 


Rearranging (4.5) gives the following expression 


V R 


R 


eff In (0. l)fl„ 2 . 3026£\ 


(4.6) 


Using (4.6), the upper 3dB frequency, / , , can now be expressed as a function 


of the rise time, t„. 


A = 5 


2.30 0.366 


"V'e ff 


2tt , 


'R 


R 


(4.7) 


Details of how the rise time is actually measured are given later in this 
chapter. 

To illustrate the effects of frequency response upon the total waveform, 
consider the ideal collimated Lyman-a signal of Figure 4.4(a), which has 
the following Fourier scries representation, 


" 2T + 


2 T A tt kW \ f 2tt kt \ 

l ~ 2 ~ 2 ~ I 1 - COS — J (cos 1 

k=l t\ k W ' T / \ T / 


(4.8) 


where T is the period of revolution of the payload and W is the width of 
the triangle pulse. In terms of the harmonic frequencies, 2 tt k/T, of the 
Fourier scries, and the rise time, the transfer function is given by. 



' f t) (ri) 


— r -1 
+ 1 


(4.9) 


which has the following magnitude and phase characteristics. 



(4.10) 


Applying these characteristics to the triangular pulse train gives the 
following infinite series representation for the output, 



(4.11) 


As an example. Figure 4.6 gives a plot of the output function, (t ) , for a 
rise time of, t B = 1.6 ms, a period of, T = 167 ms, and a pulse duration of, 

n 

U = 40 ms, as approximated on a digital computer using the FORTRAN program 
titled "FOUR".* Notice that the peak values of the triangle pulses arc 
reduced by 3.5 c o. Since an accurate peak value is needed for aspect 
information, a correction factor would be used in the final analysis for 
this example. 

4.4 Calibration and Testing 

The calibration and testing procedures for the Lyman-ot electronics are 
divided into the following three categories; (1) zero offset adjustments, 

(2) rise-time measurements and buffer-amplifier compensation, and (3) full- 
scale sensitivity measurements and selection of calibration resistors. A 
detailed description of these procedures is given in this section, preceded 
by a preliminary discussion of the rise-time measurements. 

To measure the rise time of the circuit a square wave of some pcak-to- 
peak voltage, 2 V, is applied to the input and then the duration of a 
transition from -VG.G^/Rf. to 90" of VG.G^/R is measured at the output, 

Gji.7,, is the gain of the circuit as described in Section 4.2 and R ^ is the 
output resistance of the square-wave generator. 


*A listing of the program FOUR is given in Appendix I. 
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Assuming a circuit gain of G^(J^ * 25 V/nA as in the two eclipse 
payloads and a generator output resistance of R 0 ■ 50 fi, the full-scale 
input voltage will be determined. The input resistance of the Lyman-a 
electronics, which is characterized by the inverting input of the Keithley 
302 electrometer, is given by, 


R. = R ♦ -f = 50 « ♦ 
1 o A 


2.5 x 10 0 


2.5 x 10 12 


(4.12) 


where A is the open- loop gain of the electrometer and is specified by the 
manufacturer. The current which produces a full-scale output of 5 V is 
given by, 

I ps = 5 V/G 1 J : = 5 V/(25 V/nA) = 0.2 nA (4.13) 

Therefore, to produce a full-scale swing (0 to 5 V) at the output requires 
an input square wave peak-to-pcak amplitude given by, 


7_ e = I vc R. = (0.2 x 10“ 9 A) (2.5 x 10 5 fi) = 50 uV 

ro to i 


(4.14) 


which is too small of a signal to be generated by a conventional waveform 
generator. To produce a square wave of such a small magnitude, a triangular 
wave is differentiated by coupling the output of a standard waveform 
generator to the input of the Lyman-ot circuit with a 10 pF capacitor as 
shown in Figure 4.7. The peak-to-pcak amplitude, V„, of a 100 llz triangular 
wave which produces a 50 yV peak-to-pcak square wave at the input of the 
Lyman-a electronics will now be determined. For positively sloped transition 
of the triangular wave the (displacement) current through the 10 pF capacitor 
is given by, 

t si n A l ^ rr 1 r 


I c- c Jt 


- c ft' 


(4.15) 


where f is the fundamental frequency of the triangular wave, V is the voltage 
across the capacitor, and I ^ is the capacitor current. To produce a full- 
scale output of 5.0 V requires a capacitor current of 


- r 

R. a FS 

t 


(4.16) 


The capacitor voltage, V, in terms of the triangle wave amplitude, I 
is given by, 

V - V - V 

v T FS 


(4.17) 


«< A *’ 


WAVEFORM 

GENERATOR 



LYMAN -a 


ELECTRONICS 


i R . 





TO SCOPE 
5 VOLTS P-P 
SQUARE WAVE 


18.1020 : A = 10 V p-p 

18 1021/22 : A= .1 V p-p 


Figure 4.7 Test setup for measuring rise time. 
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Combining equations (4.15), (4.16), and (4.17) gives the following expression 
for the amplitude of the triangular wave. 



(4.18) 


To complete the example, a peak-to-pcak triangle wave amplitude of, 

,-9 




(0.2 x 10 A) 


2(100 Hz) (10 x 10" 12 F) 


0.1 V 


(4.19) 


would be required to produce a full-scale swing at the output of the Lyman-a 
electronics. This is well within the capabilities of a conventional waveform 
generator. 

A description of the rise time measurement procedure using the above 
scheme for generating the small magnitude square wave is given in Table 4.3, 
Calibration Procedure 2. Also in Table 4.3, Calibration Procedure 1 gives 
the zero-offset adjustment, and Calibration Procedure 3 gives the full-scale 
sensitivity measurements and selection of calibration resistor. Figures 4.8 
and 4.9 give the circuit details of the Lyman-a electronics and the 
component layout, respectively, and should be referenced when reviewing the 
calibration procedures. Following the calibration procedures and figures 
is the calibration record, Table 4.4, for flights 18.1020, 18.1021, and 
18. 1022. 

Since each IIV deck has two identical Lyman-a circuits, one for the 
collimated detector and one for the uncollimated detector, a circuit will 
be referred to as either the collimated electronics or the uncollimated 
electronics Also, due to the extreme sensitivity of the electrometer 
amplifier to external electric fields, all measurements require the circuit 
shield, i.e., the aluminum case, to be in place and the hole in the bottom 
of the deck adjacent to the photometer mounting bracket must be covered 
with a grounded conductive material. 

4.5 Analysis of the Two Detector System 

Independent aspect information is not obtainable from conventional 
solar (or lunar) aspect sensors during the darkness of totality. However, 
it is possible to use two Lyman-a ion chambers, oriented differently with 
respect to the solar direction, to give aspect information along with 
corrected flux data. 
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Table 4.3 Calibration procedures. 

Calibration Procedure 1: Zero-Offset Adjustment 

1. Short-circuit the electrometer feedback resistor, R^, in both the 
collimated and uncollimated electronics sections. 1 2 3 

2. Connect a DVM to the collimated circuit output. 

3. Adjust the 50 kO potentiometer of the collimated electronics to produce 

0.00 V at the (collimated) output. 

4. Repeat steps 2 and 3 with the uncollimated electronics. 

Calibration Procedure 2: Rise-Time Measurements and Buffer Amplifier 

Compensation 

1. Disconnect the cathodes of the output protection diodes of both the 
collimated and uncollimated electronics to allow a negative swing at 
the outputs. 

2. Connect the test setup shown in Figures 4.7 using the uncollimated 
electronics keeping the 10 pF capacitor inside the circuit shield and 
using a shielded cable from the signal generator to the input pin. 

3. Adjust the amplitude of the input triangle wave to produce a 5 V 
peak-to-peak square wave at the output. 

4. Measure and record the rise time, which is the duration of a 0 to 90 
percent positive transition of the 5 V peak-to-peak square wave. 

5. Select the maximum value of C^^ that will not cause an increase in 
the rise time value obtained in the previous step (being sure to 

turn off the circuit when changing C.,^ to avoid damaging the Keithley 302). 

6. Repeat steps 2 through 5 with the uncollimated electronics. 

7. Reconnect the cathodes of the output protection diodes. 


Calibration Procedure 3: Full-Scale Sensitivity Determination and Selection 

of Calibration Resistor 

1. Connect the test setup shown in Figure 4.9 with the collimated electronics, 
making sure the 1.05 GO calibration resistor is inside the circuit shield 
and using a shielded cable from the precision voltage source of the 
calibration resistor. 


2. Adjust the precision voltage source to get 5 V at the collimated output 
and record the voltage. 


3. Disconnect the precision voltage source. 


4. 


5. 


6 . 


Apply a voltage equal to the value of the system voltage which is 
approximately 30 V to the +30 V CAL input of the collimated electronics 


Select a value for R^,, 
sure to turn off the^ 


to produce +5 V at the collimated output 
circuit when changing fl^^) . 


Repeat steps 1 through 5 with the uncollimated electronics. 


(being 


3 

5 



Schematic diagram of a Lyman- a electronics hoard. 
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Figure 4.9 Test setup for full-scale sensitivity measurement. 
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Table 4.4 Calibration record for flights 18.1020, 18.1821, and 18.1022 
for the collimated (C) and uncoil imatcd (U) detectors. 



18.1070 C 

18.1020 U 

18.1021 C 

18.1021 U 

18.1022 C 

t R 

(msec) 

1.2 

1.0 

1.0 

1.6 

1.3 

C 741 

(pF) 

820 

820 

820 

750 

820 

y 

/ FS 

(V) 

21.2 

20.4 

0.21S 

0.218 

0.197 

t fs 

(nA) 

20.2 

19.4 

0.208 

0.208 

0.188 
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To demonstrate this concept, consider a simple arrangement in which two 
ion chambers are oriented at 30° and 9 )° to the spin axis, as shown in 
Figure 4.10. For this preliminary analysis assume that the sensors have 
an aspect snesitivity which is represented by a cosine function.* With the 
sun at an angle o to the spin axis of the rocket the signals from the two 
detectors arc, 

I = I cos (a - 30°) 

= Iq cos (90° - a) » Iq sin a (4.20) 


where I ^ is the signal at normal incidence.** 

Solving these two equations gives the angle a and I_ in terms of the 
two detector outputs, and I 2> 


a 


tan 


( 1 5775 

VV 1 : ' °' 5 


) 


(4.21) 


and 




■ hh * J 2 


(4.22) 


The physical system also has a maximum allowable range of values for 
the angle a from which J ^ can be determined. Since I ^ and I 2 a re always 
positive the following constraints are placed on the angle a, 

I = cos (a - 30°) > 0 Z=>-60° < a < 120° 

I n = sin(a) > 0 zZ> 0 < °t $ 180° (4.23) 


thus restricting tne angle a to values within the interval [0,120°]. 

In general, any two detector aspect detection system has a range of 
usable aspect angles limited by the orientation of the sensors. Referring to 
Figure 4.10, an increase in the allowable range of aspect angles of the 
example system could be facilitated by decreasing the included angle between 
the two detectors. As an example, if detector 2 was to remain at 90° from 
the spin axis, the range of usable aspect angles would be given by the 
interval [0, 90° + a'], where angle a' is the orientation of detector 1 
and is in the interval [0, 90°]. It should be noted that a trade-off exists 


*The actual aspect sensitivity is determined from the flight data and is 
discussed in section 7.3.3. 

**J Q will vary during the flight and is proportional to the solar Lyman- a 
radiation. 
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between the range of usable angles and angular sensitivity of the system, 
i.e., to large of a range would severely decrease the resolution of this 
aspect detection system. 

t 

i Theoretically with equations (4.21) and (4.22) find data from the two 

detectors, rocket aspect and solar Lyman-a intensity can be determined, 
respectively. However the equations were derived based on the assumption 
of a cosine aspect sensitivity. The actual aspect sensitivity of the 
detectors can be closely approximated by the function (cosinc) n where n 
is somewhere in the interval 1 to 2. The actual value of n is determined 
from the flight data and is discussed in section 7.3.3. 
i The slowly varying angular response of a detector favors an accurate 

flux determination but works against an accurate determination of aspect 
angle. To create a steeper angular variation a collimator could be added 
to the detector thus making it more desirable for aspect angle measurements. 

The detectors in the eclipse payloads arc used for the determination 
of both aspect angle and Lyman-a flux. Therefore, the two detector 
configuration chosen for the eclipse payloads consists of one detector, 
without a collimator, directed toward the probable direction of the sun 
(to be used primarily for flux determination) and the other detector, with 
a cylindrical collimator, directed perpendicular to the spin axis (to be 
used primarily for flux determination) and the other detector, with a 
cylindrical collimator, directed perpendicular to the spin axis (to 
be used for angular determination and measurement of flux from the horizon). 
An explanation of how the cylindrical collimator affects the angular response 
of the detectors will follow. 

Consider a simple system consisting of a detector with a window of 
diameter 2 a fitted with a cylindrical collimator of the same diameter and 
length l, as shown in Figure 4.11. The percent area of illumination of the 
collimated detector surface as a function of the angle of incidence 8 will 
now be derived. From Figures 4.11 (a) and (b) it can be seen that equations 
(4.24) are true. 

d - Z tan 8 

a cos <p = j d (4.24) 

Combining these two equations gives the following definition for the angle <{>, 


1 


PROJECTION - 
OF APERTURE 


' d 1 
2o- 

d-X TAN 


ION CHAMBER 
WINDOW 



DETECTOR 

SURFACE 


AREA OF ILLUMINATION 


o COS <£= 1/2 d 


J AREA=o /4(2<£-SIN 2<p), 




WHERE <£= COS i /2a TAN y 
90°- 


Figure 4.11 Geometry of cylindrical collimator of diameter 2 a 
and depth 2. 
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cos 


s tan 6 


0 < * 


tr 

< T 


(4.25) 


where $ * n/2 for normal incidence. From Figure 4.11 b it can be seen that 
the area of illimination is four times the shaded area of Figure 4.11 c. 

The total area of illumination is found by evaluating the following surface 
integral , 


it/ 2 
• 

AREA OF _ . 

ILLUMINATION * 

ir/2-4> 


a 

r dr dQ 

* 

acos4>/sin6 


= a“ (2$ - sin24>) (4.26) 

Finally, dividing (4.26) by the total surface area of the detector, tta* - , 
gives the following expression for the percentage of the window surface that 
is illuminated. 


PERCENT AREA 
ILLUMINATED 


(2<p - sin 2d>)/tT 


(4.27) 


As an example. Figure 4.12 gives a plot of (4.27) for l - 2a and 
0 < 6 < 50°. Notice that the response is cut off for angles greater than 45° 
and the nearly linear variation of illuminated area with respect to angle 
of incidence. Angular sensitivity of the collimated detector can be found 
by forming the product of the uncollimated response and equation (4.27). 

The response of the two eclipse payload detectors is plotted in 
Figure 4.13 for two possible aspect angle sensitivity functions, cos 6 and 
cos‘ e. Detector 1, which is oriented at 60° to the spin axis has an 
angular variation given by, 

Ij/Iq = cos n (&j) = cos(a - 60°) (4.28) 


where 6^ is the angle of incidence for detector 1, a is the angle of the 
rocket spin axis with respect to the sun, and n = 1 or 2. Detector 2, which 
is oriented at 90° to the spin axis, has an angular variation given by. 


J 2 /J 0 


cos n (B 2 ) 


/ PERCENT AREA') 
ILLUMINATED / 


= cos n (90° - a) x (2$ - sin 2$)/* 






* 


(4.29) 


PERCENT OF AREA ILLUMINATED 


Figure 



ANGLE OF INCIDENCE (/ 3 ) 

4.12 Variation of percentage of area 

illuminated with angle of incidence, 
P, for l = 2a. 
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where 6- is the angle of incidence for detector 2, <f> * cos'* [tan(90° - a)], 

4m 

and 0 < $ < tt/2 . By examination of Figure 4.13 it can be seen that the 
latter variation is closely approximated by the function (a/45 - 1) for a 
in degrees and for 1 < n < 2. This near linear response of the detector 
is an important outcome of using the cylindrical collimator. 

For flights 18.1020, 18.1021, and 18.1022 an approximately ±5 9 « 
processional motion of the payload was expected at an average sun angle of 
a = 65°. From Figure 4.13 (for n = 2), the ratio would vary from 0.25 

to 0.51 uniquely determining the sun angle throughout the flight. The 
final arrangement of the detectors chosen for the payloads is shown in 
Figure 4.14. A photograph showing a close-up of the two Lyman-a detector 
openings as well as the visible light sensor assembly is given in Figure 4.15. 




Figure 4.14 Configuration of the two Lyman-oc detectors 

with collimated detector at 90° and uncollimated 
detector at 60°, with respect to the rocket spin 
axis. The length (Z) of the collimator is 1/2 in 
and the diameter (2a) is 3/8 in. 




Close-up view of the payload showing the collimated, and, 
below it, the uncollimated Lyman-a detectors. The visible 
light solar sensor is located above the two Lyman-a 
detectors. 
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5. SOLAR VISIBLE RADIATION 

A visible light detection system is included in each payload to provide 
an indication of when the rocket enters and leaves the shadow of totality 
(for the two eclipse flights) and to provide a check on aspect data while 
the rockets are outside of the shadow. 

5.1 Sensor and Ma8k 

The detector consists of a solar cell, sensitive to radiation in the 
visible spectrum, mounted behind a back plate with a small apcrature in 
the center which is located a fixed distance behind the face plate, as 
shown in Figure 5.1. The face plate is everywhere opaque except for an 
arrangement of a slot and four holes that can transmit visible light to the 
solar sensor. The solar cell signal is then amplified to not more than 
+5 V to maintain comparibility with the telemetry system. 

Figure 5.2 shows the face plate and corresponding signals for four 
different rocket angles. Notice that the angles 0° and ±30° are easily 
identified by recognizing a group of three pulses indicating 0°, or a group 
of two pulses indicating ±30°. The angles 0° and 30° provide a quick check 
of other on-board aspect instrumentation and can be used to derive an 
in-flight calibration curve for the detector. 

5.2 Electronic Circuit 

Figure 5.3 shows the electronics associated with the visible radiation 
detector. The solar sensor transducer, S5-30, is a photo-voltaic silicon 
device produce by Solar System, Inc. Amplification of the solar signal 
is provided by the 741C operational amplifier configured as a non-inverting 
amplifier. Selection of the resistor marked "select" in Figure 5.3 is 
discussed later. The Zener diodes at the power supply inputs and the 
resistor-diode combination at the operational amplifier output limit the 
range of signal to the 0 to 5 V of the telemetry system. 

5.3 Calibration 

The calibration procedure for the visible light sensor consists of an 
adjustment of the amplifier gain and the generation of two calibration 
curves, one giving the response for varied light intensity and the other 
giving the output for a variation in the angle of incidence of incoming 
light. 



CONVERTOR 


Figure 5.1 Exploded view of visible light detector. 
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The calibration setup requires a special attachment for mounting 
the sensor assembly to a photographic tripod, a power supply, a voltmeter, 
and a test source of visible light. The sun is used as the test source and 
therefore the calibration must be performed on a clear sunny day in a 
location with access to direct sunlight. 

The special tripod attachment provides pre-calibrated deviation in 10° 
increments of the zenith angle of the detector with respect to some reference 
zenith angle. The reference angle for the calibration procedure is that of 
the sun. To put the sensor in the reference position, the tripod is first 
adjusted so that the veritcal slot in the faceplate is perpendicular to 
earth's surface using a carpenter's level as a guide. Then the reference 
angle is set by holding a blank sheet of paper behind the sensor assembly 
and fine adjusting the vertical angle of the tripod to produce as sharp as 
possible shadow on the plank sheet. In this position the path from the sun 
to the center of the faceplate should be perpendicular to the faceplate and, 
by definition, the rocket angle in this case would be 0°. Using the adjust- 
ment on the special tripod attachment the sensor can now be tilted in 10° 
increments. Once the reference position is determined, the amplifier gain 
adjustment and calibration measurements arc performed as quickly as possible 
to avoid erroneous data due to the apparent movement of the sun in the sky. 

With the sensor in the reference position, the amplifier gain in set 
by choosing the appropriate value for the resistor marked "select" in 
Figure 5.3. For the pre-eclipse flight the resistor is chosen to give an 
on-scale output for full light at normal incidence and for the two eclipse 
flights the resistors are chosen to give approximately 4 V output for 3.2 
percent of full light at normal incidence. The following resistor values 
were chosen for the three payloads. 

18.1020 R = 820 0 

18.1021 R = 560 ft 

18.1022 R = 510 ft 

Once the gain is properly set, the first calibration measurements are 
made by placing different neutral density filters in front of the faceplate 
and then recording the corresponding voltage outputs. Using kodak Written 
neutral density filters numbers 0.5, 1.0, and 2.0, data points arc taken 
at 0.1, 0.32, 1.0, 3.2, and 10 percent transmissions. Table 5.1 gives the 
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relationship between filter numbers and degrees of transmission for the 
filter combinations used in the calibration. Once the five data points 
(six in the case of pre-eclipse flight) are taken, they arc plotted on 
semi-log paper. Figure 5.4 gives least-squares linear fits to the three 
sets of data taken for the instruments on flights 18.1020, 18.1021, and 
18.1022. 

The second set of calibration data, the angular variation, is taken 
at 1.0 percent of full sunlight and then at 0.1 percent of full sunlight. 
Readings arc taken in pairs at 0°, ±10°, ±20°, ±30°, *40°, ±50°, and the 
average of each pair is then normalized and plotted. Figures 5.5 and 5.6 
give the angular variation for the three payloads for 1.0 percent and 0.1 
percent of full sunlight, respectively. Notice that the sensor has a 
higher angular sensitivity at lower light levels. 



NORMALIZED SENSOR OUTPUT 



normalized sensor output 



0 ±10 ±20 ±30 ±40 ±50 

ANGLE OF SENSOR 


Figure 5.0 Angular variation of solar sensor for T = 0.1 %. 


66 






6. PAY l.OAU INTEGRATION 

6.1 Telemetry System 

The telemetry system of the payload is used to relay data from the 
various experiments down to the ground. The system used in the present 
rocket flights is an FM/FM system where each data signal is assigned to 
a channel using the IRIG (Inter Range Instrumentation Group) proportional 
bandwidth system. With this system sixteen channels of information are 
transmitted to the ground station on a single carrier frequency using a 
frequency division multiplexing scheme. Each data signal is fed into a 
voltage controlled oscillator with a center frequency corresponding to one 
of the IRIG channels. The voltage controlled oscillator will accept a 
signal with a 0 to 5 V swing, which produces a +7.5 percent deviation of 
the center frequency at the output. The outputs of the voltage controlled 

■ 

oscillators are mixed together and are then fed into an FM transmitter for 
transmission to the ground station. 

A block diagram of the portion of the on board telemetry system 
pertaining to experiments discussed in this report is given in Figure 6.1. 

Since the bandwidth of an individual channel varies according to its center 
frequency, the signals are transmitted on IRIG channels that have information 
bendwidths wider than the anticipated highest frequency component of the 
signal. The FM/FM system provides simple and economical data transmission 
from small sounding rockets. 

The telemetry system also encompasses the ground based instrumentation 
which is used to recieve and record the incoming telemetry signal. The 
signal is picked up with a tracking antenna and fed into a FM receiver. 

The receiver output is stored directly on magnetic tape along with a time 
code and a 100 kHz reference signal. Provisions arc also made for real-time 
inspection of the incoming data during the flight. The real-time signal is 
decoded, using phase-locked loop frequency discriminators and low-pass 
filters, and is presented on a chart record together with the time code. A 
block diagram of the ground based instrumentation is given in Figure 6.2. 

The oscilliscope and the chart recorder are used to monitor the operation 
of the experiments while the rocket is in flight. 
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Figure 6.1 


Rocket telemetry system for the UV experiment. 
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Figure 6.2 Ground telemetry system [ Evans and Smith 













69 


6.2 Control and Power Systems 

The control and power system provides power for the payload instrument- 
ation and controls certain occurances which must be excuted during a flight. 
Only the functions performed by the control and power system that pertain 
to the experiments covered in this report will be discussed. 

The payload instrumentation is powered by a 30 V rechargeable battery 
th 3 t consists of twenty 2 Ampere-hour zinc-silver oxide cells packaged in 
a 3 " x 4" x 7 ” aluminum container. Two power busses are derived from the 
battery; a 30 V line directly from the battery and a -30 V line supplied 
via a regulated dc/dc converter module. 

Provisions are made to power the payload from an external power supply 
before a launch to preserve battery charge during the pre-flight checkout 
of the payload. This external supply maintains a full charge in the battery 
and also powers a heating element in the battery compartment. The heater 
is turned on only if the payload temperature is to low for optimum battery 
performance. 

Two minutes before launch time, a relay is triggered to switch the 
payload to internal power. As a safeguard against this relay malfunction- 
ing, a network of altitude sensitive switches (baroswitches) are arranged 
to switch the payload to internal power at an altitude of 12 km, regardless 
of the status of the relay. The baroswitch network also provides a 30 V 
calibration signal to the Lyman-a electronics while the rocket is in the 
range of 40 to 79 km. The purpose of this calibration signal is discussed 
later. 

The instrumentation is protected during the actual launch of the rocket 
by two doors covering the sensor openings. Once the rocket has attained 
sufficient altitude (50 km) these doors must be released. At a predetermined 
elapsed time, a mechanical timer, which is started by the thrust of the 
rocket motor, ignites an explosive charge to jetison the doors from the 
payload. The door-ejection circuit is armed at 21 km by the baroswitch 
network. A block diagram of the pertinent features of the control and power 
system is shown in figure 6.3. 

6.3 Mechanical Layout 

The aluminum framework of the University of Illinois payload consists 
of 12 inch diameter, 1/4 inch thick decks supported by 4 rigid side rails. 
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The entire assembly is enclosed in a 3/8 inch thick cylindrical shell with 
ejcctable doors covering two opposing openings in the side of the payload. 

The doors are ejected during the flight to expose the data collecting 
instrumentation to the outside world. The location of the doors on the 
University of Illinois payload and the location of the University of Illinois 
payload with respect to other instrumentation is shown in Figure 6.4. 

Figure 6.5 gives a cutaway view of the University of Illinois payload showing 
the location of the experiments discussed in this report. Concluding this 
chapter is Figure 6.6 which is a photograph of the fully assembled University 
of Illinois payload section, taken prior to installation. 
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Figure 6.4 The general arrangement of the three 
payloads for the eclipse operation. 
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Figure b.5 Arrangement of the University of Illinois experiments 
in the payload. 




Figure 6.6 The lower section of one of the Nike Tomahawk payloads 
including the X-ray, Lynan-a and solar sensors, the 
partially extended booms carring the probes and the 
particle detectors, and supporting instrumentation. 



ORIGINAL PAGE lb 
Of POOR QUALITY 



75 


7. DATA PROCESSING 

7.1 Introduction 

The signals representing the outputs of the various experiments are 
such that they require processing before they can be presented as useful 
data. A computer is used to perform most of the processing tasks. The 
recorded telemetry signal is first demultiplexed then digitized and 
re-recorded in a digital form. The recording of the digitized data is 
provided in an IBM format from the Wallops Island facility. To take 
advantage of the high speed of the CDC 175 computer at the University of 
Illinois, the IBM formatted data are copied to yet another tape in a format 
compatible with the CDC machine, thus providing both IBM and CDC formatted 
tapes . 

This chape sr describes, in detail, the processing of the raw data and 
discusses further processing to be performed in the future. 

7.2 Telemetry Calibration 

At the beginning of each IBM -format ted data tape are five records 
representing five output levels from the FM discriminators that were used 
to decode the data from the composite frequency division multiplexed 
telemetry signal. These levels represent band center, ±3.75", and ±7.50*1. 
deviation from the center frequency of each channel. Ideally, they would 
correspord exactly to 0.0, 1.25, 2.50, 3.75, 5.00 V output from the payload 
instrumentation. In practice, this is not always the case and sometimes a 
correction must be applied to acquire the true levels. This correction is 
discussed later in this chapter. 

The calibration levels for the present experiment, as read by the 
FORTRAN program CALIB [Zirmerman and Smith, 1980], are given in Table 7.1. 

7.3 Processing of Lyman-a Data 

The recorded signals from the collimated and uncollimated photometers 
are processed to derive Lyman-a absorption profiles and ri «.'t aspect 
information. This processing is broken into the following categories, 

(1) Establishing the reference levels 

(2) Peak detection 

(3) Detector aspect calibration 






Table 7.1 Calibration levels from the digital data tapes 
as read by the program CALIB. 


Instrumentation Channel 

-7.50% 

-3.25% 

BAND CENTER 

+3.75% 

+7.50% 

18.1020 Lyman-a (U)* 

201.1 

111b. 1 

2038.7 

2957.3 

3878.1 

Lyman -a (C) 

109.7 

1073.0 

2036.3 

3001.3 

3964.8 

Solar Sensor 

145.0 

1985.3 

2030.0 

2976.5 

3925.4 

18.1021 Lyman-ct (U) 

203.4 

1118.5 

2041.4 

2960.3 

3880.9 

Lyman-a (Cl 

111.8 

1074.8 

2037.8 

3002.4 

3966.0 

Solar Sensor 

149.8 

1091.0 

2037.9 

2985.3 

3935.1 

18.1022 Lyman-a (U) 

194.7 

1113.6 

2040.7 

2963.2 

3887.5 

Lyman-a (C) 

112.2 

1075.2 

2039.7 

3005.5 

3970.8 

Solar Sensor 

152.5 

1095.1 

2041.3 

2989.4 

3940.1 


*U = Uncollimated detector 


C = Collimated detector 




7.3.1 Reference levels. Two reference voltages for each Lyman-a 
photometer are available on the digital tapes for verification of the 
calibration of the payload telemetry system. Unlike the demodulation 
operation which is performed in the laboratory with exact calibration levels, 
the modulation of the telemetry signal takes place within the payload and 
is therefore vulnerable to any variation in system performance that may be 
caused by extreme temperatures and/or mechanical shock. Thus there is a 
need to establish two reference voltage levels during the flight. These 
levels represent 0.0 and 5.0 V outputs from the Lyman-a electronics. The 
0.0 V reference is taken from the Lyman-a data just prior to the ejection 
of the doors and the 5.0 V reference is generated by the Lyman-a electronics 
beginning 23 seconds after the launch and has a 6-second duration. 

The recorded 5.0 V signal was noisy due to vibrations caused by the 
thrust of the Tomahawk stage and was not used. It is recommended that in 
future flights which use this type of instrumentation that the 5.0 V 
calibration signal be generated during a time interval of minimal mechanical 
vibration. Fortunately, the Lyman-a output before the door ejection was 
quite smooth and usable zero reference signal was obtained. The zero 
references for flights 18.1020, 18.1021, and 18.1022 are given below. 


Flight 

Number 

Detector Type 

Col limated 

Uncollimated 

18. 1020 

80.4 

155.1 

18.1021 

235.1 

00 

CM 

NO 

18.1022 

142.5 

245.9 


Once the zero references are determined, the telemetry system 
calibration levels of Table 7.1 can be normalized to the Lyman-a output. 

This is accomplished by shifting the five calibration levels the same amount 
so that the value corresponding to -7.50 percent deviation will be equal to 
the zero reference. 

As an example, consider the channel corresponding to the uncollimated 
detector of flight 18.1020. The zero reference for this channel is 155.1 
and from Table 7.1, the -7.50 percent deviation level is 201.1. Therefore, 



to normalize the calibration levels of this channel each must be reduced 
by 46 units, resulting in the following data conversion levels. 
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155.1 

1070.1 
1992.7 
2911.3 

3832.1 


0.0 Volts 
1.25 Volts 
2.50 Volts 
3.75 Volts 
5.00 Volts 


With these data conversion levels, the data points of the digital tapes 
can be translated into voltages using linear interpolation between levels: 
an FM transmission scheme inherently provides a linear output. For example 
a recorded data point of 5000 would correspond to a voltage level of 3.87 
volts. The data conversion levels for all three flights along the two 
additional quantities, the shift magnitude and the percent correction, are 
given in Table 7.2. The shift magnitude is the absolute value of the 
difference between the zero reference and the -7.5 percent deviation 
calibration level, and the percent correction is 100 times the ratio of 
the shift magnitude to the overall deviation of the calibration signal 
from -7.5 to +7.5 percent. With the aid of Table 7.2 the raw data can be 
converted to voltage levels with the percent correction indicating the 
degree of uncertainty of the data. 

7.3.2 Peak detection. The spinning motion of the payload sweeDs the 
view of the lyman-cx photometers past the sun producing a periodic output 
signal. This signal appears as a train of rounded pulses from the 
uncollimated photometers. Typical segments of Lyman-a waveforms are given 
in Figure 7.1. The maximum values of both types of pulses produced is of 
interest for determining rocket aspect and Lyman-a absorption profiles. 

These pulses will be referred to as "peaks" and their maximum values will 
be referred to as "peak value". 

The detection of peaks and the determination of the peak values is 
performed by the FORTRAN program UVPEAKS , Appendix I. It takes input from 
the CDC formatteo tapes and outputs peak values a? ng with the corresponding 
elapsed times from launch. It also determines background Lyman-a values 
which are average signal levels between peaks. This background Lyman-a 
values which are average signal levels between peaks. This background 
signal represents the scattered Lyman-a detected by the photometers when 
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the sun is not within their direct viewing angle. A detailed description 
of how the program UVPEAKS works will now be presented. 

The CDC formatted tapes are arranged in records containing one real 
number followed by an array of 5000 integers. The 5000 integers represent 
5 channels of 1000 data points each, with every set of five consecutive 
integers containing one data point from each channel. The data points of 
channel n, where 1 < n < 5, arc given by every fifth element of the 5000 
integer array starting with the n th element, i.e., channel n is made up 

of all elements with indices in the set A, where A = (n + 5j : 

je { MO [0,999] }) . The real number at the beginning of each record 
represents the elapsed time from launch, which corresponds to the first set 
of 5 data points of the record. It should be noted that the elapsed time 

from launch is actually the elapsed time from ignition of the rocket rounded 

off to the nearest second. By doing this each record can be identified by 
a whole number of seconds and if an exact time is needed, the elapsed time 
from launch can be scaled accordingly. The exact launch times as taken 
from the chart records are 16:52:02, 16:52:01, and 16:54:01 UT for flights 
18.1020, 18.1021, and 18.1022, respectively. 

The program UVPEAKS processes the Lyman-a data one record at a time. 

A specified channel of each record is copied into the last 1000 locations 
of an 1100-element array called DATA. The first time a data record is 
read in, the first 100 elements of DATA are zero. For all subsequent 
readings of data records the first 100 elements of DATA contain the last 100 
elements of the previous data record. This overlap allows detection of 
peaks whi *h span two data records. 

To detect a peak, the program scans all data points that are above a 
certain threshold that is specified by the user. The purpose of this thres- 
hold is to avoid detection of false peaks from within the background signal. 
An adaptive threshold option can also be specified by the user which allows 
the program to set the threshold equal to some fractional value of the 
magnitude of the previously detected peak. This fractional value remains 
constant throughout the run and must be entered by the user. It is some- 
times necessary to test run the program to select the proper fraction. The 
adaptive threshold option is useful for when the background signal varies 
considerably during the flight. Before discussing how a peak is actually 
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detected two definitions are necessary. An "increase" is defined as a set 
of ten consecutive data points which have at least five positive differ- 
ences between adjacent points that are less than the integer IJUMP. The 
default value for IJUMP is 100, but if necessary the user can enter a 
different value. The purpose of IJUMP is to provide noise immunity for the 
detection algorithm. Similarly a "decrease" is defined as 10 consecutive 
points with at least five negative differences between adjacent points 
which are less than the integer IJUMP. 

The program UVPEAKS will scan the data looking for an "increase". Once 
an "increase" is found the program will look for a "decrease" starting with 
the data point that is the sixth element of the "increase". If the 
"decrease" is found the program will search for and count the local maximums 
within the data points bounded by the sixth element of the "increase" and 
the sixth element of the "decrease". If one local maximum and no local 
minimums are found the local maximum is considered a peak value. If more 
than one local maximum and/or more than zero local minimums are found, the 
program will apply a smoothing algorithm until either a peak value of the 
peak is identified or until the algorithm is applied ten times. The 
smoothing algorithm is given by the following difference equation.* 

t/( n) = [x(n - 1) + .'*(n) ♦ x(n + l)]/3 

If a peak value is found, it is written to the output file along with the 
corresponding elapsed launch time to the nearest millisecond. Also, the 
number of filtering iterations used to obtain that peak value and the number 
of data points between the previous and present peak values are written to 
the output file. These numbers give an indication of the reliability of a 
specific data point. If no local maximum is found the program will proceed 
to search for the next "increase" starting with the tenth data point of the 
"decrease". The peak value identification algorithm is summarized in the 
flow diagram of Figure 7.2. 

The program UVPEAKS has been used to determine the peak values of the 
Lyman-ot data taken during flights 18.1020, 18.1021, and 18.1022. The program 
parameters used to generate these data are given in Table 7.3. The final 

*Notice the non-causal nature jf this smoothing algorithm. This is 
possible because the data for the entire flight is available from the 
digital tape in a random access manner. 


SCAN THE DATA 
UNTIL AN "INCREASE' 
IS IDENTIFIED 


CONTINUE SCAN UNTIL 
A "DECREASE" IS IDENTIFIED 


SCAN THE DATA BETWEEN 
THE "INCREASE" AND THE 
"DECREASE", INCLUSIVE, FOR 
THE CONDITION* ONE LOCAL 
MAXIMUM AND ZERO LOCAL 
MiNIMUMS. 


WRITE THE LOCAL 
MAXIMUM VALUE 
TO THE OUTPUT 
FILE 


APPLY SMOOTHING ALGORITHM 
TO THE DATA IN 
QUESTION 

TEST FOR II th SMOOTHING 
ITERATION 


Figure 7.2 Flow diagram of the peak value identification 
algorithm used in the program IJVPF.AKS . 
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Table 7.3 Parameters used for executing UVPEAKS for each detector output. 


Detector 

IJUMP 

ITHRSH 

RATIO 

IS. 1020 uncollimated 

100 

1?\. 

— 

IS. 1020 collimated 

200 

90 \ 

0 . 1 

18.1021 uncollimated 
18.1021 collimated 

200 

100 

345 ! 

/' 

235 

0.5 

0.5 

18.1022 uncollimated 

400 

400 

0.5 

18.1022 collimated 

400 

250 

0.5 
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plots are included in Chapter 8. It should be noted that these data 
are presented as a function of flight time, but could be given as a function 
of flight time, but could be given as a function of altitude using the 
flight trajectory' record to convert flight time to altitude. 

7.3.3 detector aspect calibration. The Lyman-a photometers inherently 
exhibit a variation in output level with angle of incidence of incoming 
photons. A function fit to this angular variation is derived using 
experimental data obtained during flight 18.1021. 

Before investigating the angular variation of the photometer, a 

preliminary expression which gives the included angle between a line normal 

to the surface of the photometer window and a line from the window to the 

sun in terms of the aximuthal and elevation angles of the two linear paths, 

will be derived. The included angle will be referred to as y, he aximuth 

angles will be called $ . and <j> , and the elevation angles will be identified 

s p 

as © s and 0^ , where the subscripts s and p refer to the path to the sun and 
the normal to the photometer window, respectively. A pictoral representation 
of these anglei is given in Figure 7.3. Consider the right triangle in 
Figure 7.3 with hypotenuse c. Applying the Pythagorean relationship to this 
triangle gives, 


2 

c 


2 

c “ + (b sin 6 
o v s 


a sin 6 ) 
P 


2 


(7.1) 


The law of cosines applied to the triangle defined by the three sides labeled 

a , b , and c gives, 
o o o & 


c = a + b 

- 2 a b cos(4> - <J> ) 

(7.2) 

0 0 0 

o o s p 

In terms of lines a and b 

the sides a and b are, 
0 o 


a = a cos 0 
o p 


(7.3) 

b = b cos 0 
o s 


(7.4) 


Substituting (7.3) and (7.4) into (7.2), and then substituting (7.2) into 

2 

(7.1) gives the fol’owing expression for c , 

(7.5) 

2 2 2 

c = a + b - 2ab sin 0 sin 0 - 2ab cos 0 cos 0 cos(tf> - $ ) . 

s p p s s p 


SUN 



Figure 7.3 Pictorial representation of the angles pertinent to the Lyma.i-cx 
detector aspect calibration. 
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2 

Using (7.5) to represent c , the law of cosines applied to the triangle 
defined by the sides labeled a, b, and c gives the desired expression, 


cos y = sin 0 sin 0 ♦ cos 0 cos 0 cos(4„ - 4) 

s p s p s p 


(7.6) 


which defines the angle of interest, y, in terms of the two linear paths 
a and b. 

The "peak values", as determined by the program UVPEAKS, represent the 

detector output for whenever the azimuth angle of a line orthogonal to the 

detector window is equal to the azimuth angle of a linear path from the 

rocket to the sun, i.e., whenever 4! = 4 . Since the elevation angle of the 

P s 

sun is essentially constant throughout the flight (270), a plot of this 
data would vary as a function of the photometer elevation angle, 0^, and 
the Lyman-a intensity. For determining the angular sensitivity of a 
detector, the variation due to Lyman-a intensity is eliminated by considering 
only the detector output for when the payload is in a region of the atmosphere 
wher'' no absorption of Lyman-a takes place. This region of constant Lyman-a 
flux is generally considered to be altitudes above 100 km. 

The collimated detector response is modeled by, 


I (y) = I Q cos P (y) (7.7) 

where p is some real number between 1.0 and 2.0 to be determined by function 
fitting I (y) to appropriate pulses of the detector output signal. A 
discussion of the selection of pulses suitable for modeling the detector 
angular sensitivity will follow. 

Consider the aspect position defined by 0 = b and 4 = 4 . For this 

s p s p 

oiientation of the payload it can be seen from (7.6) that y = 0 which 
indicates that the photometer window surface is orthogonal to a linear path 
from the detector to the sun. By observation of a plot of peak values, such 
as the plots in Chapter 8, it can be determined if the detector has passed 
through this orthogonal position and if it has the location where the 
payload passes through this position can be identified. 

Due to the processional motion of the payload, the peak v'aluc plots of 
Chapter 8 exhibit a consistent periodic variation throughout the flight. 

If the detector passes through the orthogonal orientation during a flight, 
an additional pair of smaller variations will be present at one extreme of 


88 


The processional cycle caused by the detector passing through the point of 
its maximum output.* The maximum values of these smaller variations occur 
when the detector passes through the orthogonal position. Therefore, to 
locate the or u'.oconal positions of a given flight, one would look for che 
maximum values of the smaller variations. 

If the orthogonal aspect portion occurs while the payload is above 
the Lyman-a absorbing region or does not occur at all, an alternate aspect 
measuring device (e.g. a spin magnetometer) must be used to indicate when 
the payload is at the closest point to the orthogonal position, while it 
is still above the absorbing region and then a scaling factor must be used 
when determinin and the angular sensitivity. Since the output of an 
uncollimated detector varies minimally for small angular deviations, a 
highly accurate alternate aspect sensor is not needed as long as the 
detector comes within a few degrees of the orthogonal aspect position. 

As an example, consider a flight where the uncollimated detector does 
not pass through the orthogonal aspect position but docs come within 2° at 
some point above the Lyman-a absorbing region, as indicated by a spin 
magnetometer. The angular sensitivity model (7.7) would be scaled down by 
a scaling factor, cos^(2°), where 1 < p < 2, before it is fitted to a pulse. 
For p = 1 the scaling factor is 0.999 and for p = 2 the scaling factor is 
also 0.999. If the spin magnetometer was inaccurate by 2° and the actual 
payload aspect was 4°, the scaling factor should have been somewhere 
between cos(4°) = 0.998 and cos“(4 <> ) = 0.995 which indicates a maximum 
scale factor error of only 0.4 percent. From this example, it can be seen 
that a high accuracy is not required of the spin magnetometer and that for 
small deviations from the orthogonal aspect position, the scaling factor is 
essentially equal to 1 and is therefore not needed. 

Once the point in the trajectory is found where the detector either 
passes through or is within a few degrees of the orthogonal position, the 
angular sensitivity model (7.7) can be fitted to the pulses to determine the 
proper value of p. An explanation of how this was done for the data of 
flight 18.1021 will follow. 

* A good example of this type of plot is the plot of peak values of the 
collimated detector output of flight 18.1020 in Chapter 8. 


**■ 
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From the spin magnetometer output, it was found that the closest the 

uncollimated detector window came to being orthogonal to a line from the 

detector to the sun was 3° at an altitude of 132 km. and 194.5 s into the 

flight. Since cos(3°) s cos 2 (3°) = 1.0 a scaling factor is not needed. 

Five pulses from around 194.5 s into the flight were chosen for the 

curve fit. A plot of these pulses is given in Figure 7.4. Notice the 

considerable background level from when the sun is not directly within 

the viewing angle of the detector and the steep sides of the pulses caused 

by the edges of the opening in the shell of the payload. 

Since the detector response is symmetrical and a smooth plot is 

desirable, each of the five pulses was divided into two parts at its axis 

of symmetry and then a point for point average of the resulting ten sets 

of data was determined. The background level was subtracted from the 

average variation and also the time scale was converted to an angular scale 

of 0° to 45°, where 0° corresponds to normal incidence. To convert from a 

time scale to an angular scale equation (7.6) was used with 0 § = 30°, 

0 = 27°, $ =0° (reference), $ = f (t) = 360° (t/r), and T = 210 ms, where t 
P P ^ 

is the spin period of the payload. Rearranging (7.6) and inserting these 
values gives the following equation for angular deviation in terms of 
elapsed time (in milliseconds) from the center of a pulse. 

y = arc cos [0.227 + 0.772 cos(£/0.583)] (7.8) 

A plot of the resulting average angular variation of the detector, denoted 
/( y) , is given in Figure 7.5 (a). 

Since the angular sensitivity is to be modeled by cos n (Y), the following 
equation must be true 

/( Y) h cos n (Y) (7.9) 

Taking the natural logarithm of both sides gives the following parametric 
equation in terms of the angle of incidence, y. 

ln[/(Y)] “ n ln[cos(Y> J (7.10) 

Since the opening in the payload shell cuts off direct Lyman-a penetration 
at 21° (see Figure 7.5), the slope of a least squares linear fit to the 
locus of points defined by, 

{ (ln[/(Y)], ln[cos(Y)]) : 0 < y < 21° }, (7.11) 
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which gives a value of n ■ 1.1414 « *^2. Figure 7.6 shows the data points in 
relation to the linear approximation and Figure 7.5 (b) gives the angular 
response model, cos' (Y)» in relation to the actual angular sonsitixity. 
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8. FLIGHT PERFORMANCE 

8.1 Eolipae Civomatanaea 

The general location of the temporary launch site established by NASA 
is indicated in Figure 8.1. The launch site is near Red Lake, Ontario and 
the exact geographic coordinates of the launchers are 50.9029 °N, 93.4583 °W. 

Two Nike Tomahawks (18.1021 and 18.1022) carrying the instrumentation 
described in this report were launched on 26 February 1979, during the 
eclipse. One other ( 1 Q . 1 0?0) . carrying similar but less sensitive solar 
radiation experiments, was launched on 24 February 1979 to obtain background 
data and to check out the launch facilities. 

The launch times for the two rockets in totality were based on a 
predicted trajectory having an effective launch elevation angle in the 
range 83 to 86 degrees, an azimuth angle of 15 degrees, and an apogee in the 
range 134 to 138 km. The actual trajectories were close to those predicted 
with apogees of 132.6 and 132.3 km for Nike Tomahawks 18.1021 and 18.1022, 
respectively. Apogee for the pre-eclipse launch (Nike Tomahawk 18.1020) 
was 130.5 km. 

The trajectories of Nike Tomahawk 18.1021 launched at 1652:00 UT and 
Nike Tomahawk 18.1022 launched at 1654:10 UT on 25 February 1979, are shown 
in Figure 8.2. The eclipse circumstances at the position of the rocket arc 
indicated along each trajectory. For times when the rocket is outside of 
totality, the number given is the percentage of the disc that is visible. 
Where the rocket is in totality the time since second contact (i.e. the 
start of totality) is given. 

It can be noted that the first of the two rockets enters totality at 
T + 60 s, at an altitude of 61 km, on ascent and exits at T + 300 s, at an 
altitude of 68.5 km, on descent. The second rocket remained in totality 
from launch until T + 106 s, at an altitude of 105 km, on ascent. 

The circumstances confirm that, as planned, one rocket was in the D 
region just after the start of totality and the other just before the end 
of totality. 

Another perspective on the eclipse circumstances is presented in Figure 
8.3. Here the apparent sun-moon distance, in solar radii, is plotted against 
time. The magnitude of the eclipse is 1.0403 (i.e. the ratio of the mean 
diameters of moon and sun) so that a sun-moon distance of 0.0403 solar radii 


RED LAKE' 



ipse of 26 February 1979 
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represents the edge of totality and 0 is center of totality. This f igui s 
gives a picture of the rocket position ir. terms of radial distance from the 
center of totality. 

A complementary view of the trajectories is shown in Figure 8.4. This 
is a polar plot of the sun-moon distance against the angular position 
relative to the sun's North point, and gives a picture of the rocket position 
viewed along the axis of the shadow. It is a remarkable coincidence that 
the relative motion of the shadow and the first rocket (18.1021) is such as 
to keep the rocket in totality for such an extended period of time (rather 
longer than was anticipated). 

1.2 Solar Radiation Experiments 

Preliminary observations of the output signals from the detectors 
indicate that good direct 3nd background Lyman-a data, as well as rocket 
aspect, will be obtainable using appropriate data processing techniques. 

The signal from t' l pre-eclipse flight exhibits a progressive decay starting 
at approximately 7 ♦ 100 s from launch. Since apogee for this flight 
occurres at T + 182 s, this decay indicates a gradual degradation of 
sensitivity during the flight, which will require additional processing to 
correct. 

At the time of the completion of this report (August 1979) the signals 
from both Lyman-a detectors of each flight have been successfully processed 
using the algorithm described in Chapter 7. Plots of the maximum values of 
the pulses as well as the background signals for both the collimated and 
uncollimated detectors for the three flights are given separately in 
Figures 8.5 to 8.7. 

For the pre-eclipse flight, Nike Tomahawk 18.1020, the background 
level in both detectors is too small to be seen in Figure 8.5. Recall that 
the sensitivity of the detector experiment of the pre-eclipse flight was 
less by a factor of 100 than that used in the eclipse flights. Since the 
detectors are exposed to the full sun (i.e. no eclipse) on that occasion 
the incident flux is constant while the rocket is above the absorbing 
region of the atmosphere. In fact the recorded signal is not constant. It 
shows a period variation, associated with the precessional motion of rocket, 
and unexpectedly, a gradual decrease in signal, indicating a degradation 
of the detectors themselves. Since this is an effect of the high 
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(unattenuated) flux the actual Lyman- a flux can be dcterrai^‘•‘' , by extrapolating 
back to the time that the rocket ex^ts the absorbing region, at an altitude 
of about 110 km. 

The rapid increase in signal in the initial part of the data represents 
the decreasing absorption of the Lyman-a flux and can be converted to a 
concentration profile of molecular oxygen. 

Neither detector on Nike Tomahawk 18.1021 show evidence of degradation, 
presumably attributable to the relatively small Lyman-a flux in totality 
(roughly 0.5"o of the flux from the uneclipsed sun). The symmetry of the 
data, in time, is consistant with the symmetrical nature of the trajectory 
of the rocket with respect to the umbra (sec Figure 8.4). The processional 
motion of the rocket is evident in the oscillating nature of the signal. 

As expected the uncollimatcd detector shows less sensitivity to payload 
orientation than does the collimated detector. The included angle of 
the precession cone for this flight, determined from the aspect magnetometer, 
is 20 degrees. 

Nike Tomahawk 18.1022 was launched during totality and exited the umbra. 
As seen in Figure 8.7 the signal increases, with a superimposed oscillation 
again due to the processional motion, until it reaches full scale (which 
represents about 1 percent of the flux from the uncclipscd sun). These 
data, together with those of 18.1021 allow the brightness of the solar 
chromosphere in Lyman-a to be determined. Also of interest in Figure 8.7 
is the variation of the intensity of the scattered Lyman-a signal. This has 
not yet been explained, even qualitatively. 
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9. CONCLUSION AND RLCOMM LNDAT IONS FOR FUTURF. WORK 

This report has described the instrumentation and preliminary results 
of a group of solar radiation experiments flown in sounding rockets before 
and during the eclipse of 26 February 1979. The over-all design of the 
structure of a scientific payload and the post-flight data processing tasks 
have also been discussed. 

Originally it was planned that results from both the X-ray and Lyman-a 
experiments would be included in this report. Due to an unusually high 
flux of energetic particles on the day of the eclipse, the hard X-ray data 
from the two eclipse launches was contaminated and is therefore not presented. 
However, it is fortunate that satellite observations on that day show that 
no significant level of hard X-rays was produced by the sun at the time of 
the eclipse. 

The Lyman-a data has been processed to some extent and is presented 
in this report. Maintaining the goals of this paper, a detailed description 
of the Lyman-a data processing tasks, already performed, has been given 
along with motivation for the final processing of the data. A summary of 
the overall Lyman-a data processing requirements will follow. 

The purpose of the Lyman-a data processing is to derive rocket aspect, 
direct lyman-a intensity, and scattered background Lyman-a intensity from 
the recorded output currents of the two Lyman-a photometers. The data 
processing is divided into the following seven tasks. 

(1) Modelling of the angular response of the uncollimated 
detector. 

(2) Modelling of the angular response of the collimated 
detector. 

(3) Determining the maximum values of the output pulses 
from the two photometers to be used in the derivation 
of both the di ect Lyman-a intensity and the rocket 
aspect profile ;. 

(4) Determining the background level between the pulses 
of both signals to be used in the derivation of both 
the direct Lyman-a intensity and the rocket aspect 
profiles. 
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(5) Deriving the rocket aspect profile from the maximum 
values obtained in (3). 

(6) Deriving the direct Lyroan-a intensity profile by using 
the rocket aspect profile to normalize the data of (3). 

(7) Derive the background Lyman-a intensity profile by 
applying a diffusion model to the data of (4). 

A functional fit has been made to the angular response of the uncolli- 
mated detector in Nike Tomahawk 18.1021. It was found that the function 


7 


/2 

cos 


(X) 


(9.1) 


is a very good approximation to the detector output, where 7 () is the output 
current of the detector at normal incidence, and 7 is the output current for 
incoming radiation at an angle of incidence X. Details of the curve fitting 
to the uncollimated detector output have been given in this report. 

The collimated detector is closely approximated by a linear function 
and requires a linear curve-fitting technique. Since the uncollimatcd 
detector does not approach a position of normally incident radiation, two 
degrees of freedom must be considered when fitting the function to an output 
pulse. 

A FORTRAN program titled UVPHAKS has been used to identify the maximum 
values of the output pulses and the background levels between the output 
pulses for both detectors of all three flights. The theoretical base and 
use of the program has been covered in this report. Outputs of the program 
have been plotted and are also presented. 

The peak values of the signals from the Lyman-u detectors are also used 
to determine the rocket aspect. Since the Lyman-a intensity varies with 
altitude, the output currents from the detector must vary as a function of 
altitude as well as aspect angle. The outputs in terms of aspect angle, X 
and altitude, are 


Ij(M) 


V*> 


cos (X) 


( 9 . 2 ) 


/ 2 (*,X) = («X + b) 


( 9 . 3 ) 


where 7 . and 7 . are the maximum values of the uncollimatcd and collimated 

1 4m 

detector outputs, repscct ively, and m and b are to be determined by curve 


fitting to the collimated output pulses. Taking the ratio of (9.2) and 
(9.3) gives the following equation. 
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cos (x) 
m x ♦ b 


(9.4) 


Since m and b are known constants and the ratio I./I. (.h,x) is easily 
calculated, the aspect angle x can be determined, as a function of altitude, 
by solving (9.4) numerically for the aspect profile, x(^)» using a digital 
computer. 

With the aspect profile available, the direct Lyman-a intensity profile 
is easily determined by solving (9.2) for (fc,x) and inserting the 
corresponding values of x(h) • The background Lyman-cx intensity profile 
is determined in a similar manner. 

As mentioned in the text, the orientation of the Lyman-a detector is 
quite critical in an aspect determination application. In future applications 
of this technique, it is advised that care be taken in positioning the 
detectors to minimize the data processing and for independence from other 
aspect measuring devices. For example, it was necessary to use the spin 
magnetometer data to indicate when the detector approached a position of 
normally incident solar radiation. If the Lyman-a detector had passed 
through this position it would have been identifiable from the Lyman-a data 
independent of the spin magnetometer output. If this was the case, the spin 
magnetometer would have served only as verification of that position as 
opposed to being the sole indicator. 

In general, the Lyman-a experiments were quite successful in that 
preliminary data indicate that accurate Lyman-a profiles, as well as rocket 
aspect, have been successfully obtained. 

As a final note, it is hoped that the results of the experiments 
discussed in this paper will aid in the understanding of the complex 
interaction between solar radiation and the various components of the 
atmosphere. 
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APPENDIX I 

Listings of the programs UVPEAKS and FOUR 
A.I.l UVPF.AKS 

PROGRAM UVPEAKS ( INPUT .OUTPUT , UVDATA, UVOUT, BOUT , TAPE1 =U VGA fA , 
1TAPE2=UVQUT, TAPES=B0U T , TAPE4= INPUT , TAPE5=0UTPUT ) 

C 

C THIS PROGRAM IS USED TO IDENTIFY AND TABULATE PEAK VALUES 

FROM UV DATA OBTAINED DURING FLIGHTS 18. 1020 , 18. 1021 , AND 18.1022. 

DATA(IIOO) = DATA RECORD WHICH IS BEING PROCESSED 

IJUMP = MAXIMUM ALLOWABLE INCREASE/DECREASE BETWEEN CONSECUTiVE 

POINTS; USED FOR NOISE ©ISCHIMIMATIOH 

ICON = NO. OF CONSECUTIVE POINTS TO BE PROCESSED DURING EACH CALL 

TO SUBROUTINES INC OR DEC 

NA = NO. OF DATA POINTS TO BE AVERAGED AT A TIME BY 

SUBROUTINE SMOOTH 

ITHRSH = INITIAL THRESHOLD FOR ADAPTIVE THRESHOLD OPTION 

THRESH = THRESHOLD CONSTANT 

RATIO = (PRESENT THRESHOLD)/( PREVIOUS PEAK VALUE) 

TRKNUM = THE TAPE TRACK TO BE PROCESSED 

T = TIME OF FIRST DATA POINT IN A RECORD 

TIP s TiME OF PREVIOUSLY IDENTIFIED PEAK 

J1 = PREVIOUS PEAK LOCATION 

INC 1 /DEC 1 = LOCATION OF FIRST DATA POINT IN AN "INCREASE" OR 

"DECHEASE" 

INCF/DECF = INDICATES IF AN " INCREASE"/"DECHEASE" WAS FOUND IN 

THE SPECIFIED SEARCH LIMIT 

NID = NUMBER OF INCREASES (OR DECREASES) USED TO .'DENTIFY AN 

"INCREASE" (OR "DECREASE") IN THE SUBROUTINE INC (OR DEC) 

SP = THE APROXIMATE SPIN PERIOD 

INTEGER DATA, DEC 1 , NID , NA , DHANL, DHANU 

INTEGER DECF , PRANL, PKANU , NA , NFi, ARHAY (5000 ) , TRKNUM , THRESH 
COMMON /ALL/ DATA(IIOO) 

COMMON /MDI/ NiD, ICON, IJUMP 
COMMON /MPDI/ THRESH 
COMMON /MPB/ T 1 P , J 1 

common /mp/ t.thryes, ratio, ithrsh 

COMMON /MB/ SP 

DATA NID, NA, THRESH, ICON, IJUMP/5, J, 425, 9, 100/ 

DATA J1 , T1P/0, 0.0/ 

DaTA STOP, ITHRSH/500. 0,0/ 

REWIND 1 
REWIND 2 
REWIND } 

INITIALIZE DATA( 1-100) TO ZERO. 

DO 10 J= 1,100 
10 DaTA(J)=0 

PRINT (5 ,*) "ENTER THE TRACK NUMBER AND THE APPROXIMATE SPiN PERIO 
ID IN MILLISECONDS;" 

READ (4 , • ) TRKNUM, SP 
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PRINT (5 , • ) "THE CURRENT PARAMETER VALUED ARE;" 

PRINT (5 , •) " THRESHOLDS" , THRESH, " , IJUMP = " , IJUMP, " , ST0P=",ST0P 
PRINT (5, # ) "DO YOU WANT TO CHANGE ANY PARAMETERS ?" 

READ 20, PAR YES 
20 FORMAT (A6) 

IF ( PAR YES . £U. "N" ) GO TO SO 
PRINT (5 , • ) "ENTER DESIRED VALUES" 

READ (4,») THRESH, IJ UMP, STOP 

print (5 f * ) "thresholds" , thresh, " , ijump=", ijump,", stop=",stop 

SO PRINT (5,») "ADAPTIVE THRESHOLD OPTION ? (Y OR N)" 

READ 40, THRYES 
40 FORMAT (A6) 

THE ADAPTIVE THRESHOLD 0P' r .0N WILL SET A NEW THRESHOLD EQuhL TO rt 

CERTAIN FRACTION (SPECIFIED BY RATIO) OF THE PREVIOUS PEAK VALUE. 

IF (THRYES. EQ."N") GO TO 50 

PRiNT (5,*) "ENTER THRESHOLD RATIO AND INITIAL THRESHOLD;" 

READ (4 , * ) RATIO, THRESH 
ITHRSHsTHRESH 

READ ONE RECORD FROM THE TAPE INTO ARRAY. 

50 READ (1) T, ARRAY 

IF (EOF(I).NE.O.O) STOP 
iF ( T . GT .STOP) STOP 

OUTPUT THE TIMES OF EVERY TENTH RECORD TO INDICATE TO 

THE USER THAT THE TAPE IS BEING READ. 

iF ( AMOD ( T , 1 0 . 0 ) . EU . 0 . 0 ) PRINT (5,*) "T= ",T 

READ THE TRACK SPECIFIED BY TRKNUM FROM ARRAY INTO 

DATA( 100-1 100) . 

DO 60 L= 1,4996, 5 
L 1 = 10 1 + ( (L-1 )/5 ) 

L2=TRKNUM+L-1 
60 DATA(L1 )=ARRAY(L2) 

INCIs-NID+1 
DECIs-NID+1 
JJ=-SP».8 

C IF A PEAK WAS FOUND SET LOWER SEARCH LIMIT OF SUBROUTINE 

C INC TO A NUMBER OF LOCATIONS PAST THE PREVIOUS PEAK 

C WriiCH IS EQUAL TO 80* OF THE APPROXIMATE SPiN PERIOD 

C IN MILLISECONDS. IF A PEAK WAS NOT FOUND SET THE LOWER 

C oEARCH LIMIT TO THE TENTH ELEMENT OF THE PREVIOUSLY 

C DETECTED "DECREASE" . 

C 

fO IRANL=D£C1+N1D+(SP«.8)+JJ 
JJ=0 
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IF ( IRANL.GT . 1090) GO TO 80 


SET UPPER SEARCH LIMIT FOR SUBROUTINE INC TO LOCATION 1100. 

IHANU= 1 1 00 

CALL INC ( IRANL , IRANU , INCF, INC 1 ) 

IF (INCF.EQ.O) GO TO 80 
IF (IIOO-(INCUNID).LT.IO) GO TO 80 

SET LOWER SEARCH LIMIT FOR SUBROUTINE DEC TO INCUNID. 


DHANL=INCUNID 

SET UPPER SEARCH LIMIT FOR SUBROUTINE DEC TO LOCATION 1100. 

DHANU= 1 1 00 

CALL DEC ( DRANL , DRANU , DECF , DEC 1 ) 

IF (DECF.EQ.O) GO TO 80 

SET SEARCH LIMITS FOR SUBROUTINE PEAK. 

PRAN. xINCUNID 
PRANU=DECK10-NID 

IF ((PRANL.LT. 2). OR. (PRANU.GT. 1099)) GO TO eo 
call peak (pranl,pranu,na,nfi) 

GO TO to 
80 CONTINUE 

c j 

C MOVE LAST HUNDRED ELEMENTS OF DATA(IIOO) INTO DATA( 1-100). 

C I 

DO 90 K= 1,100 
K 1 = 1000+K 

90 DATA(K)=DATA(K1 ) 

GO TO 50 
END 

C 

SUBROUTINE PEAK ( PRANL, PHANU , NA , NFI ) 

C J 

C THIS SUBROUTINE WILL SEARCH FOR A ''PEAK" WITHIN THE LOCATIONS OF 

C DATA(IIOO) SPECIFIED BY THE SEARCH LIMITS PRANU & PRANL. IF 

C NECc-SSAHY , SUBROUTINE PEAK WILL CALL THE SMOOTHING SUBROUTINE 

C SMOOTH UP TO NINE TIMES TO IDENTIFY A "PEAK". IF A "PEAK" IS 

C FOUND, THE PEAK VALUE IS WRITTEN TO THE OUTPUT FILE UVOUT ALONG 

C WITH THE CORRESPONDING TIME AND NUMBER OF SMOOTHING ITERATIONS. 

C THEN A CALL IS MADE TO SUBROUTINE BACK AND THE RESULTING 

0 BACKROUND VALUE AND CORRESPONDING TIME ARE WRITTEN TO OUTPUT 

C FILE BOUT. IF NO "PEAK" IS IDENTIFIED AFTER NINE CALLS TO SMOOTH, 

C CONTROL IS RETURNED TO THE MAIN PROGRAM AND NOTHING IS WRITTEN TO 

C THE OUTPUT FILES. 

C PRANL = LOWER SEARCH LIMIT 

C PRANU = UPPER SEARCH LIMIT 

C NFI = NUMBER OF FILTER ITERATIONS 




l 
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C NOB = NUMBER OF OAT* POINTS BETWEEN PREVIOUS PEAK ANO 

C PRESENT PEAK 

C NLM = NUMBER OF LOCAL MAXIMUMS 

C J1 = PREVIOUS PEAK LOCATION 

C J2 = PRESENT PEAK LOCATION 

C T1 = TIME CORRESPONDING TO PRESENT PEAK VALUE 

C TIP = TIME COfiESPONDING TO PREVIOUS PEAK VALUE 

C JB = NUMBER 0/ LOCAL MAXIMUMS 

C 

INTEGER DATA, NFI, NLM , PRANL , PHANU , NLM , NDB , THRESH 
COMMON /ALL/ OATA(IIOO) 

COMMON /MPDI/ THRESH 
COMMON /MPB/ TIP , J 1 
COMMON /MP/ T.THRYES, RATIO, ITHRSH 
COMMON /PB/ NDB.T1 
NFI=0 
10 NLM=0 
JB=0 
J 3=0 

DO 30 1=PRANL , PRANU 

C 

TEST FOR THE MAXIMUM ALLOWABLE NUMBER OF CALLS TO 

SUBROUTINE SMOOTH 

IF (NFi.Ey.10) GO TO 140 

COUNT NUMBER OF LOCAL MXNiMUMS. 

IF ( ( DATA( 1 ) . LT . DATA( 1+1)) .AND . ( DATA ( I ) .LT. DA TA( 1-1 ) ) ) JB=JB+1 

TEST FOR CANDIDATES FOR LOCAL MAXIMUMS. 

IF ( ( DATA( I ) . GE. DATA( I+1 ) ) .AND. ( OATA( j. ) .GE . DATA( i-1 ) ) ) GO TO 20 
GO TO 30 

SET J2 EQUAL TO THE LOCATION OF THE PRESENT CANDIDATE FOR A 

LOCAL MAXIMUM. J3 IS THE LOCATION OF ThE PREVIOUS CANDIDATE FOR 

LOCAL MAXIMUM. 

20 J 2=1 


IF TWO CONSECUTIVE CANDIDATES ARE FOUND, CONSIDER THE FIRST ONE 
TO BE A LOCAL MAXIMUM AND IGNORE THE SECOND ONE. OTHERWISE THE 
LOCATION J2 IS CONSIDERED TO INDICATE A LOCAL MAXIMUM 
AND THE LOCAL MAXIMUM COUNT (NLM) IS INCREASED BX ONE. 

IF ( ( J2-J i ) . EQ. 1 ) GO TO 30 
NLM=NLM+1 

set j 3 equal to the location of the present local maximum 

and JUMP TO ThE TOP OF THE DO LOOP TO TEST THE NEXT POINT 
'WITHIN THE SEARCH LiMlTS. 
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C 

J j*J2 

30 CONTINUE 

RETURN TO MAIN PROGRAM IF NO LOCAL MAXIMUM IS FOUNO. 

IF (NLM.EU.O) 00 TO 140 

IF EXACTLY ONE LOCAL MAXIMUM AND NO LOCAL MINIMUMS ARE FOUND 

00 TO 60. OTHERWISE CALL SUBROUTINE SMOOTH. 

It ( ( NLM. £Q. 1 ) . AND . ( JB. EQ. 0) ) GO TO 60 

FOR THE FIRST CALL TO SUBROUTINE SMOOTH INCREASE THE UPPER 

SEARCH LIMIT BI &LEVEN LOCATIONS. FOR SUBSEQUENT CALLS THE 

UPPER SEARCH LIMIT IS NOT CHaNuED AGAIN. 

IF (NFI.EQ.O) GO TO 40 
C.0 TO 50 

40 PRANU=PHANU+1 1 

IF ( ( ( PHANL-40 ) . LT . 1 ) .OR . ( ( PRANU+40 ) .GT . 1100)) GO TO 140 

APPLI SMOOTHING TO A WINDOW WHICH IS 40 LOCATIONS WIDER ON 

EACH SIDE THAN THE SEARCH LIMIT. 

50 CALL SMOOTH (NA , NFI, PRANL-40, PRANU+29) 

NOW THAT THE DATA HAS BEEN SMOOTHED JUMP TO 10 TO SEARCH AGAIN 

FOR A SINGLE "PEAK". 

GO TO 10 

60 IF (J2.LE.100) GO TO 70 
GO TO 80 

CALCULATE THE TIME CORRESPONDING TO THE SINGLE LOCAL 

. .. MAXIMUM WHICH WILL BE CONSIDERED TO BE A "PEAK". 

70 T1=T-1 .0+((899.0+FLOAT(J2))/1000.0) 

GO TO 90 

80 T1=T+( ( FLOAT (J2)-101 )/1000) 

IF TWO PEAKS ARE IDENTIFIED WITHIN 0.025 SECONDS OF EACH OTHER, 

THE SECOND ONE IS IGNORED AND CONTROL IS RETURNED TO THE MAIN 

PROGRAM. 

90 IF ( ABS(T 1P-T 1 ) . LT . 0. 025 ) RETURN 

IF THE PRESENT PEAK IS THE FIRST ONE FOUND IN THE PRESENT 

RECORD GO TO 100. 


C 


IF (J1.EU.0) GO TO 100 
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C APPLY AN ADAPTIVE THRESHOLD IF IT IS GREATER 

C THAN OR EQUAL TO THE INITIAL THRESHOLD, OTHERWISE APPLY THE 

C.... INITIAL THRESHOLD. 

C 

IF ( THRYES . eaJ. " Y" ) T.4RESH=IFIX(RATI0»FL0AT(DATA(J2) ) ) 
if ( THRESH. LT.ITHRSH) TriHESHsiTHRSH 

C 

C CALCULATE THE NUMBER OF DATA POINTS FROM TrtE PREVIOUS "PEAK - ’ TO 

C THE PRESENT "PEAK". 

C 

NDB si FIX(T1 # 1000.0)-IFIX(T1P* 1000.0) 

CALL BACK (J2.BVAL.TBACK) 

WRITE (2,120) T1 ,DATA(J2) , NFi.NDB 

GO TO 110 
C 

C IF THE ADAPTIVE THRESHOLD OPTION WAS CHOSEN AND THE PRESENT 

C PEAK VALUE IS MORE THAN 501 GREATER THAN THE INITIAL 

C THRESHOLD RETURN CONTROL TO THE MAIN PROGRAM. 

C 

100 IF ((DAfA(J2).GT.(IFIX(1. 5*FLQAT(iTHRSH) ) ) ) . AND. (THR YES. EQ. 

1 "Y ")) RETURN 
WRITE (2,1j}0) f 1,DATA(J2) ,NFI 
110 J 1 =J 2 
TIPsTI 

120 FORMAT (1H ,F7. 3, 3X,I10, 3X,I2, 3X, 16) 

1 SO FORMAT (Iri ,F7. 3»3X, 110, 3X,I2) 

140 RETURN 
END 
C 

SUBROUTINE DEC (DRANL, DHaNU , DECF , DEC1 ) 

C 

C THIS SUBROUTINE WILL SEARCH FOR A ‘'DECREASE" WITHIN THE LOCATIONS 

C OF DATA(IIOO) SPECIFIED BY THE SEARCH LIMITS DnANL & DRANU. IF 

C A "DECREASE" IS FOUND SUBROUTINE DEC WILL RETURN WITH DECF*1 AND 

C WITH DEC 1 CONTAINING THE LOCATION OF THE FIRST ELEMENT 

C OF THE "DECREASE". 

C DRANL s LOWER SEARCH LIMIT 

C DRANU s UPPER SEARCH LIMIT 

C DECF =0 IF NO DECREASE IS FOUND iN SPECIFIED RANGE 

C DECF s 1 IF DECREASE IS FOUND IN SPECIFIED RANGE 

C DEC 1 = LOCATION OF FIRST DATA POINT IN A "DECREASE" 

C 

INTEGER DATA, NID, DECF, DRANL, DRANU, DEC1 , THRESH 
COMMON /ALL/ DATA (1100) 

COMMON /MDI/ NID, ICON, IJUMP 
COMMON /MPDI/ THRESH 
DECFsO 
IDATAsO 
C 

C TEST FOR AT LEAST NID DECREASES, WHICH ARE LESS THAN IJUMP, 

C IN EACH SET OF TEN CONSECUTIVE POINTS WITHIN SPECIFIED RANGE. 

C 
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10 JYES=0 

DO 40 JTEN= 1 , ICON 

IF (DATA(DHANL+IDATA+JTEN-1 ) .LT. THRESH) GO TO 40 

IF ( DATA ( DRANL+IDATA-hJ TEN ) . LT . DATA ( DRANL+IDATA+J TEN-1 ) ) GO TO 20 

GO TO 40 

IF ( ( DATA (DRANL+IDATA+J TEN- 1 )-DATA( DHANL+IDATA+JTEN ) ) . LT . I JUMP ) 

1 GO TO 30 
GO TO 40 
30 JYES=JYES*1 
40 CONTINUE 

IF (JYES.GE.N1D) GO TO 50 
C 

C THY NEXT SET OF TEN DATA POINTS ONLY IF WITHIN SPECIFIED RANGE. 

C 

IDATA=IDATA+ 1 

IF ( ( DRANL+IDATA+9 ) .GT . DHANU ) GO TO 60 
GO TO 10 

C 

C INDICATE A DECREASE WAS FOUND AND DETERMINE THE LOCATION OF THE 

c — first data point in the "decrease". 

c 

50 DECF= 1 

DEC 1 =DRANL+IDATA 
60 RETURN 
END 

C 

SUBROUTINE INC ( IRANL , IRANU , INCF, INC 1 ) 

C 

C THIS SUBROUTINE WILL SEARCH FOR AN "INCREASE" WITHIN THE 

C LOCATIONS OF DATA(IIOO) SPECIFIED BY THE SEARCH LIMITS IRANL A 

C IRANU. IF AN "INCREASE" IS FOUND, SUBROUTINE INC WILL RETURN 

C WITH INCFsl AND WITH INC1 CONTAINING THE LOCATION OF THE FIRST 

C ELEMENT OF THE "INCREASE". 

C IRANL = LOWER SEARCH LIMIT 

C IRANU = UPPER LIMIT OF RANGE 

C INCF=0 IF NO INCREASE IS FOUND IN SPECIFIED RANGE 

C INCF= 1 IF INCREASE IS FOUND IN SPECIFIED RANGE 

C INC 1 = LOCATION CF FIRST DATA POINT IN AN "INCREASE" 

C 

INTEGER DATA, NID, THRESH 
COMMON /ALL/ DATA( 11001 
COMMON /MDI/ NID, ICON, iJUMP 
COMMON /MPDi/ THRESH 
INCF=0 
IDATArO 
C 

C TEST FOR AT LEAST NID INCREASES, WHICH ARE LESS THAN IJUMP, 

C IN EACH SET OF CONSECUTIVE POINTS WITHIN SPECIFIED RANGE. 

C 

10 J YES=0 

DO 40 J TEN= 1 , ICON 

IF (DATA(1RANL+IDATA+JTEN-1 ) .LT. THRESH) GO TO 40 
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IF (DATA( IRANL+lDATA+JTEN ) .GT . DATA( IRANL+iDATA+JTEN-1 ) ) GO TO 20 

GO TO 40 

20 IF ( ( DATA( IRANl+IDATA+JTEN )-DATA( IRANL+IDATA-hJTEN-I ) ) . LT . IJ UMP) 

1 GO TO jO 
GO TO 40 
30 JXES=JX£S+1 
40 CONTINUE 

IF (J XES.GE.NiD) GO TO 50 
C 

THX NeXT SET OF TEN DATA POINTS ONLX IF WITHIN SPECIFIED RANGE. 

IDATAdDATA-fl 

IF ( ( IRANL+iDATA+9 ) .GT . IRaNU ) GO TO 60 

go ro io 

INDICATE AN INCREASE WAS FOUND AND DETERMINE THE LOCATION OF TttE 

First data point in the " increase". 

c 

50 INCFsl 

INC1=IHANL+IDATA 
60 RETURN 
End 
c 

SUBROUTINE SMOOTH (NA , NFI , PRaNL, PRANU) 

C 

C THIS SUBROUTINE WILL SMOOTH THE DATA POINTS IN DATA(IIOO) 

C THAT ARE wiTHIN THE LOCATIONS BOUNDED BX PRANL & PRANU. THE 

C CURRENT VALUE OF NFI WILL bE INCREMENTED BX 1 EACH TiME 

C THIS SUBROUTINE IS CALLED. 

C NA = NUMBER POINTS AVERAGED 

C NFI = NUMBER OF FILTER ITERATIONS 

C 

INTEGER NA , NFI , PRANL , PRANU , I , NAT , SUM , DATA 
COMMON /ALL/ DATA(IIOO) 

NF1=NFI+1 

IF (NFI.E0.10) GO TO 30 
DO 20 I=PRANL, PRANU 
NAT=NA/2 
SUM=0 

DO 10 J = 1 , NA 

10 SUM=SUM*DATa(I-NAT-kJ-1 ) 

20 DATA( I ) =SUM/NA 
30 RETUHN 
END 
C 

SUBROUTINE BACK (J2,BVAL,TBACK) 

U 

C THIS SUBROUTINE WILL DETERMINE AND WRITE A BACKRUUND VALUE 

C TO OUTPUT File BOUT EXCEPT IF TttE PRESENT "PEAK" AND THE PrlEV i. JUS 

C "PEAK" ARE NOT IN THE SAME RECORD OR IF THE DISTANCE FROM THE 

C PREVIOUS "PEAK" TO THE PRESENT "PEAK" IS NOT WITHIN A J.S 

C PERCENT DEVIATION FROM THE APROXIMATE SPlN PERIOD. 
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J1 = PREVIOUS PEAK LOCATION 
J2 = PRESENT PEAK LOCATION 
BVAL = BACKROUND LEVEL 

TBACK = CORRESPONDING TIME OF BACKROUND LEVEL 

INTEGER DATA, BVAL 
COMMON /ALL/ DATA(IIOO) 

COMMON /PB/ NDB,T1 
COMMON /MPB/ TIP , J 1 
COMMON /MB/ SP 

TEST TO SEE IF THE PRESENT AND PREVIOUS "PEAKS" ARE NOT IN 
THE SAME RECORD OR IF THE DISTANCE BETWEEN THEM IS NOT 
WITHIN 9.6 PERCENT OF THE APROXIMATE SPIN PERIOD. 

MiN=lFIX(SP-(SP # . 096 ) ) 

MAX=IFIX(SP+(SP a .096) ) 

IF ((NDB.LT.M1N).0R.(NDB.GT.MAX).0R.(J1.GT.J2)) RETURN 

CALCULATE THE TIME CORRESPONDING TO THE PRESENT BACKROUND VALUE. 

TBACK= ( ( T 1 -T IP ) /2 . ) +T 1 
KMID=J UNDB/2 
JCOUNTsO 
ISUM=0 


DETERMINE THE PHESENT BACKROUND LEVEL. 

I1rKMlD-25 
i2=KMlD+25 
DO 10 1=11,12 

ISUM=ISUM+DATA( 1 ) 

BVAL=ISUM/51 

WRITE (3,20) TBACK, BVAL 

FORMAT ( 1 H , F7. 3, 3X , 1 10 ) 

RETURN 

END 

END 
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A. 1.2 

Sample 

output of the 

program UVI’f.AKS 

1: lapsed time 

Peak 

Number 

of 


Number of data 

from launch 

value 

filter 

Herat ions 

between peaks 

107. Jib 

812 


l 


207 

10’ .568 

821 


1 


212 

107.779 

79b 


-> 


211 

tO". 990 

798 


2 


211 

108.201 

799 


n 


211 

108.409 

798 


2 


208 

108.021 

801 


« 


212 

108.832 

760 


4 


211 

109.042 

801 


1 


210 

109.253 

797 


3 


211 

109.4b3 

796 


2 


210 

109.072 

815 


i 


209 

109.883 

793 


•j 


2il 

110 093 

812 


i 


210 

110.305 

811 




212 

110.51b 

802 


•> 


211 

110.72b 

820 




210 

110.93b 

825 


0 


210 

111.147 

824 


2 

4b 


211 

111.358 

817 


2 


211 

111. 5b7 

838 


i 


209 

111.779 

847 


i 


212 

111.989 

844 


i 


210 

112.198 

871 


i 


209 

112.410 

8b0 


i 


212 

112.619 

872 


i 


209 

112.832 

889 


i 


213 

113.043 

874 


2 


211 

113.251 

89o 


2 


208 

113.404 

904 


i 


213 

113. b72 

91b 


i 


208 

113.885 

930 


i 


213 

114.095 

93b 


i 


210 

114.305 

942 


i 


210 

114.515 

955 


i 


210 

114.724 

967 


i 


209 

114.935 

983 


i 


211 

115.147 

987 


i 


212 

115.358 

997 


i 


211 

115.507 

1014 


i 


209 

115.777 

1011 


i 


210 

1 15.990 

1042 




213 

110.199 

1033 




209 

1 10.410 

1073 


i 


211 

110.021 

1058 




211 

lio.831 

1 0 ’4 


2 


210 

117.042 

1098 


1 


211 

1 17.253 

1093 


1 


211 

117.401 

1115 




208 

ir.b’M 

1128 


l 


213 

11 ’.883 

1139 


l 


209 

118. 094 

1 1'9 


0 


211 

118.305 

1160 


l 


211 

118.51b 

1180 


l 


21 1 

118.725 

1 180 




209 

118. 937 

1190 


"9 


212 

119.148 

1222 


1 


211 

119.358 

1215 


1 


210 

119. Sh9 

1265 


0 


211 

119.780 

12b4 


1 


211 

119.990 

1250 


1 


210 

1 20 . 20 0 

127b 


1 


210 
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A. I. 3 FOUR 


PROGRAM FOUR ( INPUT , OUTPUT , FILE , TAP£1=F1LE,MIK£ , TAP£2=MIKE) 

REAL FN , F , AN 
REWIND 1 

PRINT*, "ENTER NUMBER OF ITERATIONS AND RISE TIME" 

READ* ,N, TR 
N=N+1 

DO 10 IT=1 ,501 

T= . 001 *FLOAT( IT ) — . 25 1 

FN= . 12 

F=0 . 0 

DO 20 1 = 1, N 

F=F+FN 

3=FL0AT(I) 

PI = i. 1415926 
W=6.*3 

AN=(8.i5/(Pl*PI*i*3))*(1-COS(.2i95*PI*S)) 

IF(IT.EQ.I) WR1 TE(2 , 60 ) W,AN 
60 F0RMAT(F12.5,F12.5) 

THE TA = AT AN ( ( TR »S ) / . 06 1 2 ) 

20 FN=AN*COS((S*r*2.»Pl-ThETA)/.l67)/SORT(((PI*S*TR)*»2.)»27.051+1 . ) 
10 WR1TF( 1 , 30) T,F 
}0 FORMAT ( 1X,F7.i,jX,F12.5) 

STOP 

END 


